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Abstract. The coherent mixing of the current states in the superconducting weak 
link subject to a Josephson phase difference (f> and subject to an external transport 
current in the banks is one of the aims of this work. At <f> = it the nonlocal 
mixing of current states produces two vortices close to the point-contact between 
superconducting bulks. The effect of point-contact reflection in an impenetrable 
interface and effect of temperature on the vortices have been studied. It is obtained 
that increasing the reflection of the point-contact destroys the vortices while increasing 
the temperature restore these vortices. The vortex state is a new version of the 
interference between the macroscopic states and quantum tunnelling. Also, the weak 
link between unitary triplet superconductors which have /—wave and p + h— wave 
pairing symmetry has been studied from the spin and charge current-phase relation 
point of view. The main result in the second part of this thesis, is the polarization 
of the spin transport when a junction between triplet superconductors is used. It is 
observed that the spin current is the result of the misorientation between the gap 
vectors of two supercondcutors. In addition, the weak link between two bipolar non- 
unitary triplet superconductors is studied mathematically. The current-phase relations 
obtained in third part of this thesis are totally different from the junctions between 
the unitary spin-triplet superconductors and between the spin-singlet superconductors. 
The current phase diagrams which have been obtained in this work can be used to 
distinguish the symmetry of the order parameter in the crystals. 
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1. Introduction 

Superconductor is a substance that below a certain temperature (critical temperature 
T c ) conducts charge current without resistance and also repels the magnetic field. 
Superconductivity was discovered by Kammerlingh Onnes in a mercury wire at 
temperatures less than T = 4.2K in 1911 pQ. In 1957, Bardeen, Cooper and 
Schrieffer (BCS), using the concept of Cooper pairs, proposed a microscopic theory 
for superconductivity 0. Cooper had shown earlier that the ground state of a metal, at 
sufficient low temperatures, will be stable when its electrons are collected as the pairs 
[Sj. BCS theory suggests that Cooper pairs at low temperatures condense into the same 
quantum state defined by a macroscopic wave function and can travel together without 
dissipation. Also, the pair amplitude of this wave function as an order parameter in 
the Ginzburg-Landau (GL) theory of superconductivity, defines the transition between 
normal and superconducting states [4 . In the limit of T — » T c , the GL order parameter 
, is directly proportional to the BCS energy gap, A, which is the amount of energy 
for breaking a Cooper pair. This energy gap carries the information about the pairing 
symmetry of the electrons. 

There are different types of superconductivity. According to the spin states of two 
electrons in the Cooper pair, superconductor is either spin-singlet or spin-triplet. Since 
electrons are fermion, the total wave function of a cooper pair should be antisymmetric 
with respect to electrons: ^cooper-pair = r 2V( s i; s 2) where r, s are spatial 

coordinates and spins, respectively. For the antisymmetric spin-singlet state, the orbital 
part of the wave function should be symmetric under the interchange of the electrons 
and the orbital momentum should be even (I = 0, s— wave, I = 2, d— wave and I = 4, 
g— wave). 

The superconductors with simplest spin-singlet state, I = 0, are known as conventional 
superconductors. All the other type of superconductors including the spin-singlet and 
spin-triplet with I ^ 0, are unconventional superconductors. A s— wave superconductor, 
has an isotropic order parameter A(k) = Ao in the momentum space, where k is a 
unit vector pointed on the Fermi surface. Earlier superconductors which had been 
found in the elements like mercury, are conventional superconductors. Unconventional 
superconducting compounds which have anisotropic order parameter can be defined 
using the relation A(k) = 0, with summation over the Fermi surface. For the 
symmetric spin-triplet state, the orbital part of the wave function is antisymmetric and 
the orbital momentum takes the odd numbers (/ = 1, p— wave, / = 3, /—wave and 
I = 5, h— wave). Here, the terms s— wave, p— wave and etc have been used from the 
terminology of the Hydrogen atom. The pairing state of a spin-triplet superconductor 
in the spin space is represented by a three dimensional vector d(k), called gap vector. 
The gap vector determines the order parameter matrix, A(k) = i(d(k) ■ a)a y , in which 
ajs are Pauli matrices. 

Three complex components of gap vector (di, g?2, d^) over the Fermi surface are 
corresponding to three possible spin directions as follow: \spin) = (di + zd 2 )| JJ. 
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) + (— di + id 2 ) \ TT) + ^3(1 Tl) + I IT))- There are many experimental and theoretical 
works in the case of spin-triplet superconductors from which some important properties 
of spin-triplet superconductors are listed below: 

1) The spin-triplet superconductors are generally low T c superconductors, as compared 
with the high T c superconductivity in d— wave superconductors [SI El El • 

2) The relation between critical temperature and energy gap at zero temperature is 
different from that of the BCS relation for a s— wave superconductor A(T = 0) > 1.76T C 

PET 

3) In structures with spin-triplet superconductivity a spin supercurrent can flow 
generally, while in the case of the spin-singlet superconductivity, the spin current may 
exist only in the proximity system of the superconductor and a ferromagnet. This means 
that in the spin-triplet superconductors not only the charge but also the spin of electrons 
can become superffuid jH3 [TTJ H2j . 

4) Ferromagnetic superconductivity is another interesting phenomenon of spin-triplet 
state. It had been observed that all conventional superconductors were non- 
magnetic materials and it was concluded that superconductivity and ferromagnetism 
are incompatible phases. While, ferromagnetic superconductivity has been observed 
recently, in some of the triplet superconductors like: ZrZn2, UGe2 and URhGe2 

oiei. 

5) Nonunitary gap vector is another fingerprint of the spin-triplet superconductivity. In 
the nonunitary spin-triplet state, Cooper pairs may carry a finite averaged intrinsic spin 
momentum. This nonunitary state is a candidate for the B— phase of superconductivity 
in the UPt 3 compound (FigJHJ). This phase has been observed at the low temperatures 
and low magnetic fields [T5| IT6] . 

6) In addition to the temperature and magnetic field, pressure influences the phase 
transition of triplet superconductors, particularly for the ferromagnetic superconductors 

[THEE!. 

7) Another property of some of triplet superconductors in the similarity of their 
structures with some of the important high T c superconductors. For example, 
the Sr 2 RuC>4 spin-triplet superconductor is isostructural to the spin-singlet high- 
temperature superconductor LaBaCuO. Here, Sr and Ru atoms are counterpart 
of La(Ba) and Cu atoms, respectively. Consequently, investigation of this triplet 
superconductor helps to understand the singlet case [5 . 

We have to introduce some of unconventional superconducting compounds, particularly 
those which will be investigated in the following chapters of this work. The uranium 
compound UPt 3 will be studied in chapter (jU) and (JHJ). Then, Sr 2 Ru0 4: and PrOs^Sb^ 
compounds will be investigated in chapters (JU) and (JHJ), respectively. 

The first discovered triplet order parameter, p— wave pairing symmetry, has been 
observed in the Helium superffuid. Also, p— wave order parameter has been considered 
as a candidate for the superconducting state in Sr 2 RuOi by some authors [THl l2Hj . 
A famous form of p— wave pairing symmetry in momentum space is d(T, k) = 
A(T)(k x + ik y )z in which, z is a unit vector jHj. Another important category of triplet 
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superconductors is /—wave superconductivity, proposed for the pairing symmetry in 
heavy- fermion compound UPt^ in [2B 122] and the compound Sr 2 RuO A in 123] . 
The different phases of /—wave superconductivity have different order parameter 
symmetries in momentum space. For instance, the /—wave axial state symmetry 



a nonunitary /—wave gap vector in the B— phase of superconductivity (low temperature 
and low magnetic field) in UPt$ compound has been considered. For this nonunitary 



has been proposed in the momentum space in ^H]- A recently proposed spin-triplet 
state is the"(p + h)— wave" pairing symmetry which, has been considered for the 
superconductivity in PrOs^Sbn compound in [Hj. Of course, two different phases of 
"(p + h)— wave" have been observed for superconductivity in the PrOs^Sb^ |23 EE]- 
Here, A— phase is high magnetic field, high temperature phase but B— phase is low 
field, low temperature phase [2Z]- The first model to explain the properties of the 
A-phase of PrOs^Sb^ is: d(T, k) = A (T)(k x + ik y )^(l — k x — ky — kj)z, where, 
z is a unit vector. The function A = A (T) describes the dependence of the gap 
vector d on the temperature TjH]. The second model to describe the gap vector 
of the 5-phase of PrOs A Sb l2 is: d(T,k) = A (T){k x + ik y )(l - k*)z$. One of 
the most important types of unconventional superconductivity has been observed in 
the high T c superconductor in the ceramic of copper oxide by Bednorz and Miiller 
This spin-singlet superconductor has a d— wave pairing symmetry, I — 2, as 
A(k) = A(T)(A; X 2 - k y 2 ) or A(k) = A(T)2k x k y , in momentum space. This discovery 
was a real revolution in the field of superconductivity. 

Conventional and unconventional superconductors are usually differentiated by 
several properties such as anisotropicity of the order parameter in the momentum space, 
heat capacity and heat conductance, density of states, junction behavior and particularly 
nodes in the momentum space. Node is a direction in the momentum space where no 
order parameter and superconductivity is effective on the scattering electrons. 

One of the most interesting concepts in the field of superconductivity is, 
superconducting weak link. The weak link experiments are classified in the categories 
of S — I — S, S — N — S and S — c — S, in which two superconducting bulks have 
been separated by a nonsuperconducting interface or point-contact. Here, S, I, N 
and c denote superconductor, insulator, normal metal and point-contact respectively. 
The weakness of the link means that the superconducting order parameters are the 
same as the disconnected massive superconducting bulks. The first configuration had 
been investigated by Josephson, in which a thin insulator is located between the two 
superconducting bulks. In the second case, a normal layer has been sandwiched by two 
superconducting bulks. The third experiment is devoted to the geometry consisting 
two superconducting bulks which are separated by an impenetrable interface (strong 
insulator) in which a contact has been prepared for the mixing of two superconducting 
states. Josephson effect, S — I — S, has been investigated in 1962 by Josephson 




bipolar state, a gap vector 
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PU] . He predicted an electrical current flowing between two superconductors while 
they are separated by an insulator. This insulator layer is very thin and electrons 
can pass through it, even with the energies less than height of potential barrier of 
insulator. The flow of current between the superconductors in the absence of an applied 
voltage is called a Josephson current (there is a phase difference instead of the applied 
voltage) and the motion of electrons across the barrier is called Josephson tunneling. 
This tunneling phenomenon, Josephson effect, can be used in the electronic devices 
such as superconducting quantum interference device (SQUID), for detecting the very 
small magnetic flux (even a flux quantum). The Josephson effect can be influenced by 
the external magnetic fields through the external phase difference. Consequently, the 
Josephson junction can be applied to measure the extremely weak magnetic fields, in 
SQUID. The current which flows, (quantum mechanically tunnels through the potential 
barrier of the interface, has a form j(0) = j c sin0 in which, is the macroscopic 
phase difference between two superconductors and the critical current j c depends on the 
geometry of system (junction). The Josephson effect is a weak link of two massive banks 
of superconductor, Si and £2 with different phases (pi and 2 which are separated by an 
insulator. The system allows electron to exchange between the two sides of the interface 
and then establishes the phase coherence in the system as a hole. The Josephson junction 
can be considered as the mixer between the superconducting quantum macroscopic 
states. The result of mixing is the supercurrent which flows from one of the banks to 
the other and it depends on the phase difference <p = <p2 ~ <Pi across the weak link. The 
function j(0) depends on the geometry of the system. In the simplest case which had 
been studied by Josephson, the current has a sinusoidal dependence on the phase but 
there are many different kinds of the weak links in which the current phase diagrams are 
not sinusoidal. One of the most famous problems in this category is problem of Kulik- 
Omelyanchouk (S — c — 5* junction) which has been investigated in [3U]. They have 
used the Eilenberger equation [HI] and using the Green function they have obtained the 
current-phase diagrams analytically. The current phase relation depends not only on 
the manner of coupling but also on the properties of the supercondcuting bulks. Two 
coupled unconventional supercondcuting massive bulks have totaly different current- 
phase relation from the conventional superconductivity which was studied by Josephson 
and later by Kulik and Omelyanchouk. For instance D—c—D weak link which is a special 
type of S — c — S system, in which D denotes the d— wave and high T c superconducting 
bulks and c is the contact, has some new results. The spontaneous current parallel to the 
junction interface, mid-gap states resulting from the sign change of the order parameter, 
and the changing period of current j(0) to j(20) are different and new results of weak 
link between unconventional superconducting bulks. In addition, it is observed that 
the Josephson junction depends not only on the external Josephson phase difference 
resulting form the external magnetic flux but also on the misorientation angle between 
two superconducting crystals. In this work, because of similarity between our problem 
and problem of Kulik-Omelyanchouk, we use the generalized form of their formalism. 
This thesis consists of three parts. At first, we are to investigate coherent current 
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states in the superconducting weak link which has been subjected to the Josephson 
phase difference and subjected to the external transport supercurrent state (parallel 
to the junction interface) in the banks. Earlier it had been observed that at <fi close to 
7r the mixing of current states produces two vortices in the vicinity of the point-contact 
between superconducting bulks |32]- In this part we study the effect of transparency 
coefficient (potential barrier) of the point-contact in an impenetrable interface and effect 
of temperature on the vortices obtained in the paper 

The second part of this dissertation is devoted to the weak link between unitary 
triplet superconductors. These superconductors have /—wave and p + h— wave pairing 
symmetry. The former has been proposed for UPt% and Sr^RuO^ compounds and 
the later case is considered in PrOs^Sbyi complex. The spin-current in the junction 
between these unitary triplet superconductors is an important part of this work. The 
interesting case which is observed in the second part of the thesis, is the polarization 
of the spin transport using the junction between unitary triplet superconducting bulks. 
The third and last part of the thesis discusses the case of weak link between non-unitary 
triplet superconductors. The idea behind the first part of the thesis is suitability of 
this structure (vortex-like currents) for the investigation of the quantum macroscopic 
phenomena. For the real system the finite transparency coefficient (finite reflection) 
which we have investigated is more suitable than ideal transparent point-contact which 
had been considered in The second part of the thesis discusses the case of the 

junction between unitary triplet /—wave and p + /i— wave superconductors. First of all 
we know that the triplet superconductors also the high T c superconductors are created 
by a different mechanism of pairing than the phonon-electron interaction. Secondly, 
the molecular structure of the Ruthenate compound Sr 2 RuOi is the same as that of 
the high T c Cupperate superconductors which are important superconductors and we 
are interested to understand their properties. Also, the spin polarized transport has 
a valuable motivation for physicist in the field of spintornics because, the sensitivity 
and accuracy of the polarized spin-current systems can be used in the measurement 
technology. The second part of this work can be used to develop the theory of spin- 
polarized transport systems. The spin transport in the absence of the charge transport 
is an interesting case to investigate. The third part of the work has been devoted to 
non-unitary weak links which has intrinsic spin and angular momentum of systems. 

The method for investigation of these weak link experiments is the quasiclassical 
method of Green function. This method had been used by Kulik and Omelyanchouk in 
paper [3U|. They have studied junction through the point-contact between two static 
conventional superconducting bulks, but, here we have applied this formalism for the 
case of conventional superconducting bulks with the external transport current in the 
banks. Also, we have generalized the Kulik-Omelyanchouk formalism for the case of 
triplet supercondcuting bulks. We have calculated the analytical Green's function and 
then we have used that to obtain the current density. The current-phase diagrams are 
plotted and in some cases the two dimensional profile of the current has been plotted 
in the space. 
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Arrangement of the rest of this thesis is as follows. In Chapter (J5J), we review of 
the concepts which we use in the rest of the thesis. The quasiclassical approach which is 
widely used in the field of solid state physics and specially for the case of superconducting 
systems, will be reviewed. The Kulik-Omelyanchouk problem, that is the ballistic point 
contact between two superconducting massive bulks will be studied and their results 
will be reproduced. The effect of transparency coefficient for a point contact with finite 
reflection in the impenetrable interface between two bulks ( related to our problem in 
the chapter (J3J)) will be investigated and some analytical relations for this system will 
be obtained. The Josephson junction between two unitary superconductors will also be 
studied in chapter (J2J). We will generalize this approach in the Chapters (jlj), © and 
©■ 

A Josephson effect in the ballistic point contact with transport current on the banks, 
taking into account the reflection of electrons from the contact, will be investigated in 
Chapter The contact is subject to the phase difference and the transport current 
tangential to the boundary of the contact. As it was shown in [32], in the contact with 
direct conductivity at <fi = tc and near the orifice the tangential current flows in the 
opposite direction to the transport current, and there are two vortices. It is found that 
by decreasing the transparency, the vortex-like current will be destroyed. On the other 
hand, as the temperature is increased the vortices are restored. They continue to exist 
for transparencies as low as D = | in the limit of T — > T c . This anomalous temperature 
behavior of the vortices is an interesting result which have been obtained. 
In Chapter (J3J), we have studied the spin and charge current in the ballistic Josephson 
junction in the model of an ideal transparent interface between two misorientated /- 
wave superconductors subjected to a phase difference 0. Our analysis has shown that the 
misorientation and different models of the gap vectors influence the spin current. The 
misorientation changes strongly the critical values of both the spin current and charge 
current. It has been shown that the spin current is the result of the misorientation 
between the gap vectors. Furthermore, it is observed that the different models of the 
gap vectors and geometries can be applied to the polarization of the spin transport. 
In addition, it is observed that in certain values of the phase difference <p, the charge- 
current vanishes while the spin-current flows, despite the fact that although the carriers 
of spin and charge are the same (electrons). 

A stationary Josephson junction as a weak link between PrOs^Sb^ triplet 
superconductors will be investigated in Chapter Recently, the "(p + h)— wave" 
form of pairing symmetry has been proposed for the superconductivity in PrOs^Sb^ 
compound [25J. The quasiclassical Eilenberger equations are analytically solved for 
this system. The spin and charge current-phase diagrams are plotted and the effect of 
misorientation between crystals on the spin current, and spontaneous and Josephson 
currents is studied. It is found that such experimental investigations of the current- 
phase diagrams can be used to test the pairing symmetry in the above-mentioned 
superconductors. Also, it is shown that this apparatus can be applied as a polarizer for 
the spin current. 



Charge and Spin Transport in Superconducting Weak Links 



8 



In Chapter a stationary Josephson effect in a weak link between misorientated 
nonunitary triplet superconductors is studied. The non-self-consistent quasiclassical 
Eilenberger equation for this system has been solved analytically and the current- 
phase diagrams are plotted for the junction between two nonunitary bipolar /—wave 
superconducting banks. A spontaneous current parallel to the interface between 
superconductors has been observed. Also, the effect of misorientation between crystals 
on the Josephson and spontaneous currents is studied. Such experimental investigations 
of the current-phase diagrams can be used to test the pairing symmetry in the above- 
mentioned superconductors. In Chapter (J7J), the thesis will be finished with some 
conclusions. 
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2. Superconducting weak links 

When two superconducting bulks are connected to each other by an insulator from which 
the electrons can tunnel, we have a weak link. The weakness of the link means that the 
superconducting order parameters have their value in the bulks and they are the same 
as the disconnected massive superconducting reservoirs. The Josephson effect arises in 
a weak link of two separated (by insulator) superconducting bulks with different phases. 
The electrons can be exchanged between two superconducting bulks and the system (two 
bulks and contact) tends to be the phase coherent. Mixing the superconducting states 
through the contact or link causes the supercurrent from one of the banks to the other 
bank. The current is present because of the phase difference between the bulks and 
this phenomenon can be observed in the absence of any voltage. The phase difference 
between the bulks which plays a central role in the weak link phenomena can be the 
result of the external magnetic field which is surrounded by the junction and bulks. The 
Josephson phase is a kind of the Aharonov-Bohm phase [SSI- This phase is related to 
the magnetic flux which flows from the system as = ^ § A • dl, where, $o — ff is the 
quantum of flux and A is the vector potential. The current in the junction which is a 
supercurrent and is called Josephson current can be calculated from jj = j^^t, in which, 
the E is the energy of the junction. In the next sections we review some of these weak 
link experiments analytically. The quasiclassical Eilenberger equations have been used 
to investigate these weak link systems in this thesis. This method will be explained in 
Sec l2.ll In Sec 12. 21 as an application of the quasiclassical Eilenberger equation, we solve 
the problem of a conventional superconducting bulk without any contact. These results 
are exactly the same of results of the standard BCS formalism which can be exerted 
directly on the uniform bulk system. Sec l2.2l is devoted to the Kulik-Omelyanchouk 
problem and at the end of this section we have generalized this method to a system of 
a contact with finite transparency, which is the Zaitsev problem jSlj. This problem is 
as same as our problem in Chapter (3). In Sec 12.31 we have a review in the case of the 
junction between unitary /—wave superconductors which had been done in paper [Hj. 

2.1. Quasiclassical approach 

The normal metals and superconductors can be investigated using the Green functions 
|3*K] . It has been shown by Eilenberger that Gorkov equations for the Green function can 
be transformed to transport-like equations for a quasiclassical Green function [31 J . These 
are called Eilenberger equations. Two conditions for applicability of the quasiclassical 
approach are that the characteristic length scales should be much larger than the Fermi 
wavelength and energies must be much lower than the Fermi energy e-p = Tp, hereafter 
fee = 1 and h — 1 for simplicity. The Green functions which will be used in our work, 
are Matsubara Green functions written in Nambu space. They are 4x4 matrices in 
a direct-product space of particle-hole and spin spaces. The general energy integrated 
Green function in k-space is of the form g(k,r,e m ). Here e m = nT(2m + 1) are the 
discrete Matsubara energies m = 0,1,2.... The odd integer value of (2m + 1) is the 
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result of Fermion behavior of electrons jSHl- Finally, the Eilenberger equation a ballistic 
case (there is no scattering) is as follows: 



v F • Vg + 



e m a 3 + iA,g = 0, (1) 



where, v> is the Fermi velocity and 03 = 03 <8> I in which bj (J = 1, 2, 3) are Pauli 
matrices. This is a first-order differential equation for the Matsubara propagator 
g(k,r,e m ) along classical trajectories of quasiparticles. The Eilenberger equation is 
not enough to make the solution unique and so, a separate normalization condition has 
to be introduced [3T|. With a suitable choice of condition satisfied by physical solutions 
of Eilenberger equation, normalization is written as gg — 1. To give a closed system of 
Eilenberger equations and normalization conditions it should be supplemented by some 
self-consistency equations for the self-energy A which will be introduced later in suitable 
forms. Finally, the Matsubara propagator g which satisfies the Eilenberger equation, 
normalization condition, continuity across the interfaces and self-consistency condition 
can be written in the form |37j : 

« = [ 9i + Si • o- (92 + g2 ■ or) ib 2 \ 

\icr 2 (#3 + g3 • <t) idr 2 (-g 4 + g4 • cr)ia 2 J ' 

where, the matrix structure of the off-diagonal self energy A in the Nambu space is 
X _ ( ( A + d • \ 

A ~ ^ ia 2 (A* + d* ■ a) J' [6) 

The A(k) = A(— k) refers to the spin-singlet but the d(k) = — d(— k) has been 
considered for the case of spin-triplet superconductivity. Fundamentally, the gap (order 
parameter) has to be determined numerically from the self-consistency equation, while 
in some cases, we use a non-self-consistent model for the gap which is much more suitable 
for an analytical calculation. The solution of Eq. (0) allows us to calculate the current 
densities. The expression for current is: 

j (r) = AmeTN (0) ^ (v f9i (y F , r, e m )) VF (4) 

m>0 

where, (-..)^ stands for averaging over the directions of an electron momentum on the 
Fermi surface and N (0) is the electron density of states at the Fermi level of energy. 
For the case of the s— wave superconductors our Green matrix changes to a matrix as 
follows: 




9 = ~ , 5 

Y 10*93 -9i J 

and for the case of the order parameter we have: 
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The selfconsistent equation for the spin-singlet case is 

A(r,T) = 27rAT^^ 2 (v F ,r)) VF (7) 

m>0 

where A is the electron-phonon constant of interaction and (••■) v is the averaging over 
directions of yp. The other important physical variable which can be derived form this 
Green function is the density of states. It is related to the diagonal term of the Green 
matrix as follows: 

N(E) = N(0)Re[ 9l (e m -> -iE + 0)] (8) 

There are many applications of quasiclassical method and in this chapter some of them 
will be reviewed. 



2.2. Superconducting bulks 



The simplest case of superconducting system which can be investigated by Eilenberger 
equation is a single bulk of superconductors without any contact or interaction with 
other world. In this system all of superconductivity properties are uniform and spatially 
constant. Because of the uniform properties the gradient term in Eilenberger equation 
(JTJ) is zero, 



e m a 3 + iA,g 







and the bulk solutions are: 



91 



V £ m 2 + [A] 



and 



£2(23) 



^ + |A| 2 ' 



(9) 



(10) 



respectively. In conclusion of Eq. 

N(E) L 



, the density of states for \E\ ^ |A| is as follows: 

V^W (12) 

while for \E\ ^ |A| is N(E) = 0. This expression has been obtained before from the 
BCS theory directly and using the particle-hole analysis of a bulk of superconductor. 
Also, the self-consistent equation for this system is: 

A(T) 



A(T) = 2tvXTJ2 



(13) 



m>0 v-m +|A(T)|' 

which is the same to the BCS self consistent equation. This latter had been obtained 
from the second quantization and quantum field theory method. This simple self- 
consistent equation only can be solved numerically, but near the T = and T = T c 
it has been solved analytically. Close to the zero temperature the gap function varies 
in terms of temperature as: 



A(T) = A(0) - v/27rTA(0)exp-(A(0)/T) 
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and near the T = T c the gap dependence on the temperature is as follows: 

A(T)=(^) 5 v^(7^T). 

Also, for the case of the gap function at the T = we have A(0) ~ 1.76T C . So, the 
current density for this system can be calculated as follows: 

(v F ) = j (r) = AmeTN (0) ]T (v F > = 0. (14) 

m>0 + 1^1 

As it is clear from the above expression, the current density for the homogenous 
superconducting bulks is zero. 

2.3. Kulik-Omelyanchouk problem: superconducting weak link through a point contact 

We consider the Josephson S — c — S weak link as a microbridge between thin 
superconducting films of thickness 2a (look at FigCJ). The length L and width d of 
the microbridge, are assumed to be less than the coherence length £o- O n the other 
hand, we assume that L and 2a are much larger than the Fermi wavelength and 
use the quasiclassical approach. The point-contact is an ideal transparent area for the 
electrons and there is not any reflection for the electron. We choose the z-axis along the 
interface and the y-axis perpendicular to the boundary; y — is the boundary plane 
(FigJTJ. If the film thickness d <C £o then in the main approximation in terms of the 
parameter g?/£ the superconducting current depends on the coordinates in the plane 
of p = (z,y). The open form of the Eilenberger equation for the case of spin-singlet 
superconducting systems is as follows: 

V ~lH^ + iA * n92 ( n ) ~~ iAn 9Hn) = °l ( 15 ) 

(9(7 2 ( n ) 

V—Qf- + 2e m92(n) - ^n9l( n ) = ^ (16) 

- 2e m g Hn) + 2iA* n g Hn) = 0; (17) 

where, t = y/\v y \ on the Fermi surface, r] = sgn(v y ) and n — 1,2 label the left and right 
hand superconducting bulks, respectively. Using the quasiclassical approximation, we 
select the solution the for this problem as follows: 

+ a n exp (-2sQt) ; (18) 
+ 6 n exp(-2sfit) ; (19) 
+ d n exp{-2sQt) ; (20) 







9l(n) — 


~Q ' 




iA n 


92 ( n ) = 


ft 






93( n ) — 


ft 



where, s = sgn(y) and ft = ^/e m 2 + |A| 2 . By substituting in the Eilenberger equation 
(JTJ), we obtain: 

SHn) = ~S + an exp ( _2sfit ) ; ( 21 ) 
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Figure 1. Model of an ideal transparent point-contact as an orifice in the thin 
impenetrable insulating partition. 

92( „ ) = ^ + a„( i -^)exp(-2^ f ); (22) 

*« = ^ + °"(i^) eXp( - 2sSii)i (23) 

In the main approximation on the small parameter a/£o -C 1, the self-consistency can 
been ignored and the model, in which the order parameter is constant in the two half- 
spaces 

A(r,T) = A(T)exp(^) 

in which <fi is the phase difference between superconductors, can be used. Solutions of 
Eqs. ((TJ) should satisfy the continuity of solutions across the contact y = 0,\z\ < a 
and specular reflection condition for y — 0, \z\ > a. In addition, far from the contact, 
solutions should coincide with the bulk solutions. Consequently, we find the diagonal 
term of Green functions which will be used in calculation of the current density, as 
follows [3E]: 
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Because of integration over the fermi surface (J3J) the antisymmetric and imaginary part 

of Green function remains. It is as follows: 

, . . . 77IAI 2 sin f cos f , N 

Im ( 9l (y = = - A ~ L —i h 25 

e m 2 + |A| 2 cosf 



and also 



,T V (v Fffl ) = - V ' A l 2si ^ , = tanh I (26) 

^ £m 2 + I A| 2 COS t 8 V 2T 



m>0 m>0 fc m "T |^| WJ& 2 

where we have done the angular integration and we have substituted (vfv) = \ 
and tanh ^ = J2 m >o x 2 +(2m+i) 2 m ^ ne above-mentioned relation. Consequently the 
current density has an expression as follows [3UJ EH : 

/ (0) = ^SeA(T)N (0) ^ sin ± tanh f A( ^M (27) 

where, in the case of the low temperatures T — > we have / (0) = ^A- sin | (20] while 
for T -»• T c we obtain 7(0) = 2^sin0 [H], in which, R^ 1 = ±e 2 SN(0)v F is the 
Sharvin resistance of the junction in the normal state [IHj and S is the effective square 
of the contact. This means that near the critical temperature current-phase relation is 
sinusoidal like Josephson prediction, while in the low temperatures the current-phase 
relation is non-sinusoidal and we have some unusual jumps at = ir. Now, we want 
to investigate the effect of transparency coefficient and reflection of point-contact on 
the current density which has been studied before by Zaitsev [34 . Using the Zaitsev 
quasiclassical boundary conditions for the case of Matsubara Green function we have 
obtained: 

/ / u LtolAPsinf cos | , s 

Im si (y = = ' 2 2 -2-, 28 

e m 2 + |A| 2 (l-Dsin| ) 

where D is transparency coefficient which is related to the interface physical properties 
as the potential barrier against our tunneling phenomenon. The transparency coefficient 
usually depends on the direction of scattering electron velocity, but here for simplicity 
we use the constant transparency coefficient formalism. The current density exactly at 
the contact is as follows: 



7T sincbD / A(T)Wl- J D(sinf) 2 \ 

J(0) = jSeA(T)N (0) v F — Sm ^ tanh OL 2> .(29) 

4 v/l-^smf) 2 V / 

Obviously, for the low value of transparency coefficient which is called tunneling limit, 
D — > 0, the current has linear dependence on the transparency and it is a sinusoidal 
function of phase and we have I (0) = jSeA(T) N (0) vpD sin tanh (^p- j ■ So, at the 
tunneling limit and T — > 0, the current is: I (0) = jSeA(0)N (0) v F D sin (j). These 
two former expressions for current are linear functions of transparency coefficient and 
sinusoidal functions of phase 0. Also in the limit of T — > T c , the current has the form 
of : I (4>) = ^fr-SeN (0) Vp sin 0D. This means that at high temperatures the current 
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is linear function of transparency and has a sinusoidal dependence on the phase. The 
nonsinusoidal dependence of phase will be happened for both the low temperature and 
high value of transparency coefficient. Consequently, nonlinearity and nonsinusoidal 
current-phase relation are coupled with each other. For the density of states using the 
quasiclassical formalism we obtain: 



EVE 2 - A 2 

N(E) = N(0)Re[ gi (e m - -iE + 0)] = ^ - r . (30) 

E 2 — A 2 (l — Dsin | ) 

In the limit of D — > 1 we obtain the density of states of the problem of Kulik- 
Omelyanchouk, N(E) = N(0) EVE ''~ -^i, and in the limit of D -> we have a 

E 2 — A 2 cos ^ 

disconnected system and the density of states for the system tends to the density of 

1 Ve?-a 2 - 



states of the bulk, N(E) = N(0) /r? $_ A2 , as we expect from the BCS theory [3*3] . 



2.4- Charge transport in the weak link between unitary and triplet super condcutors 

In this section we theoretically study the stationary Josephson effect in a small ballistic 
junction between two spin-triplet superconducting bulks with different orientations of 
the crystallographic axes. We consider a model of a ballistic point contact as an 
orifice with a diameter 2a in an impenetrable for electrons partition between two 
superconducting half spaces. We assume that the thickness of interface, d, is much larger 
than the Fermi wavelength and use the quasiclassical approach. In order to calculate the 
charge current in point contact we use Eilenberger equations (JTJ and the normalization 
condition gg = 1. For the case of the unitary and pure spin-triplet superconductivity, 
the Green function matrix (J2J) is written in the form [SJ EJ7] : 

g=( 9l + Sl ' & * g r™ 2 V (31) 
y tcr 2 g 3 -a g A - a 2 g4 ■ aa 2 J 

Matrix structure of the off-diagonal self energy A in Nambu space is 

^AttHJ-.^T 2 )' (32) 

Also, we have: 

AM H2 % i ). (33) 

Below we consider a unitary states, for which d x d* = 0. Solutions of Eq. (JTJ must 
satisfy the conditions for Green functions and vector d in the bulks of superconductors 
far from the orifice: 

^(Too) = ; ==; (34) 

e 2 m + |di |J 



l l,2 



d(Tx) =d L2 (k)cxp( T ^). (3.-)) 
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Figure 2. Scheme of a ideal transparent point-contact in an impenetrable flat interface 
between two superconducting bulks, which are misorientated as much as a. 



where <fi is the external phase difference. Eq. has to be supplemented by the 
continuity of solutions at the contact plane (\y\ < a) and conditions of reflection at the 
interface between superconductors, remainder part of the interface (|y| > a). Below we 
assume that this interface is smooth and electron scattering is negligible. The solution 
of Eq. (0) should be used to calculate the current density. We consider a simple 
model of the constant order parameter up to the surface. We assume that the order 
parameter does not depend on the coordinates and in each half-space equals to its 
value far from the point contact which is called superconducting massive bulk. For this 
non-self-consistent model the current-phase dependence of a Josephson junction can be 
calculated analytically. In a ballistic case, the system of 13 equations for functions gi 
and gi can be decomposed on independent blocks of equations. The set of equations 
which enables us to find the Green functions are: 

77-^r + HS2d -g 3 d) 
<9g- 

V ~df + ^ d X Ss + d * xg2 ) 

V^rr + £ m g2 - igid - d x g 
dt 

V-QjT - £ rng3 + iQld ~ d X g 



0; 


(36) 


0; 


(37) 


0; 


(38) 


0; 


(39) 
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where g_ = gl ~ g4 , t = y/\v y \ on the Fermi surface and r\ = sgn(v y ). The Eqs. (|36| ) -()39| ) 
can be solved by integrating over the ballistic trajectories of electron in the right and 
left half-spaces. The general solution satisfying the boundary conditions at infinity is: 

g (n) = £zn + « n exp (-2sQ n t) ; (40) 

gi n) = C n exp(-2 S fi n t); (41) 
g („) = ™nd n + d ra x C n ^ ^_ 2sQnt) . (42) 

i L n STj\L n 

g ( n) = i<± + ia n <%- d; x c n exp ( _ 2s ^ } . (43) 

where t is the time of the flight along the trajectory, sgn (t) = sgn (y) = s, r\ — sgn (v. 



V) 



and O n = + |d„| . Index n numbers left (n = 1) and right (n = 2) half-spaces. By 
matching the solutions (j40M3j) at the orifice plane (y = 0), we find constants a n and 
C n . Exactly at the orifice plane we obtain: 



Cj = C 2 ; (45) 

zdi ^ iaxdi + di x Ci _ id 2 _ ia 2 d 2 + d 2 x C 2 _ 

Qi r]Q n + e m tt 2 V^2 -£ m 

id\ zaid* — d* x Ci id* 2 ia 2 d 2 — d 2 x C 2 



(47) 

fii ryfii — e m fi 2 ?7fi 2 + e m 

Consequently, the function g\ (0) which determines the current density at the contact is 
as follows: 

n [di • di (77O2 + £ m ) 2 - d 2 ■ d 2 (r/fii - e m ) 2 ] 

gi (0) = — 2 

di • di (77^2 + e m ) + d 2 • d 2 (77^1 - e m ) + 2d x • d 2 (77^1 - £ m ) (^2 + 

Using the gi (0) one can calculate the current density at the orifice plane j(0): 

00 „ 

j(0) = 4ireN(0)v F Tj2 / d 3 kk gi (0). (49) 

m=0 ^ 

Misorientation of the crystals produces a spontaneous current along the interface |4*H l4*2"] 
generally, as can be calculated by projecting vector j at the corresponding direction. To 
illustrate the results obtained by computing the formula (}4*%j) . we can plot the current- 
phase diagrams for the different models of the pairing symmetry and for two different 
geometries. These geometries are corresponding to the different orientations of the 
crystals in the right and left sides of the interface (FigEJ). In the right hand side of 
the interface, the ab— plane has been rotated around the c— axis and the c— axis has 
been rotated around the b— axis by a in geometries (i) and (ii), respectively. Also for 
the further calculations we need to a certain model of the gap vector d which is called 
order parameter vector. There are three models which have been successful to explain 
properties of the three phases of triplet superconductivity in UPt^ compound (FigOJ). 
For the high-temperature and low field phase, A-phase, of superconductivity in UPt 3 
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Figure 3. Superconducting phase diagram of heavy fermion UPt^ compound, in the 
magnetic field and temperature space. 



the order parameter (gap vector) has an equatorial line node and two longitudinal line 
nodes [TJlEilE]]- This state which is known as the polar state, is as follows: 

d =A zk z (hi - fcj) . (50) 

The gap vector dependence in momentum space for the low-temperature and low field 
phase, B-phase, or the axial state is as follows P IT^ ITT)]: 

d =A zk z (k x + iky) 2 . (51) 

Here, the longitudinal line nodes are closed and there is a pair of point nodes. The 
coordinate axes x, y and z are chosen along the crystallographic axes a, b and c as the 
left of FigHJ The function A (T) describes the dependence of the order parameter d 
on the temperature T. Other candidate to describe the orbital states implying the weak 
effective spin-orbital coupling in the c— phase of superconductivity in U Pt 3 compound, 
is the unitary planar state jHJ 113 HE] : 

d =A k z (x (k 2 x - fcj) + y2k x k y ). (52) 

Using these forms of the triplet order parameters, we can plot the Josephson current- 
phase relation jj{4>) = j y (y = 0) calculated from EqHHl for a particular value of 
misorientation angle a under the rotation of ab-plane, the geometry (i), and rotation 
around the normal axis y or geometry (ii). For simplicity we use the spherical model of 
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the Fermi surface. The different gap vectors have different current-phase diagrams and 
such a different behavior can be a criterion for distinguish between the different phases. 
In some cases, the Josephson current formally does not equal to zero at = 0. This state 
corresponds to a spontaneous phase difference, which depends on the misorientation 
angle a. The tangential components of current, x and z, as the functions of are not 
zero when the Josephson current is zero. This spontaneous tangential current is due to 
the specific "proximity effect" similar to spontaneous current in contacts between "d- 
wave" superconductors [I2J ESI • The total current is determined by the Green function, 
which depends on the order parameters in both superconductors. As a result of this, for 
nonzero misorientation angles the current parallel to the surface can be generated. It 
can be shown that the current-phase relations are totally different for different models of 
the gap vector. Because the order parameter phase depends on the momentum direction 
on the Fermi surface, the misorientation of the superconductors leads to spontaneous 
phase difference that corresponds to the zero Josephson current and to the minimum 
of the weak link energy. This phase difference depends on the misorientation angle 
and can possess any values. It is observed that, in the "/—wave" superconductors the 
spontaneous current can be generated in a direction parallel to the plane of contact. 
Generally speaking this current is not equal to zero in the absence of the Josephson 
current. Finally, study of current-phase diagrams of Josephson junction for different 
misorientations can be used to distinguish or demonstrate the different phases and 
different triplet gap vectors of superconductivity UPt 3 . 



Charge and Spin Transport in Superconducting Weak Links 



20 



3. Coherent mixing of Josephson and transport supercurrents 

3.1. Introduction 

The investigations of Josephson effect manifestations in different systems are continuing 
due to it's importance both for basic science and industry. A point contact between two 
massive superconductors (S-c-S junction) is one of the possible Josephson weak links. A 
microscopic theory of the stationary Josephson effect in ballistic point contacts between 
conventional superconductors was developed in jHOj- Later, this theory was generalized 
for a pinhole model in 3 He ^J03], for point contacts between "ci-wave" [%2*1 146j . and 
triplet superconductors [Hj. The Josephson effect is the phase sensitive instrument for 
the analysis of an order parameter in novel (unconventional) superconductors, where 
current-phase dependencies Ij ((f)) may differ essentially from those in conventional 
superconductors EH]- in some cases the model with total transparency of the 
point contact does not quite adequately correspond to the experiment, and the electron 
reflection should be taken into account. The influence of electron reflection on the 
Josephson current in ballistic point contacts was first considered by Zaitsev [33]. He 
had shown that reflection from the contact not only changes the critical value of 
current, but also the current-phase dependence Ij (<fi) ~ sin (0/2) at low temperature 
which has been predicted in [30J. The current-phase dependence for small values of 
transparency, D >C 1, is transformed to the Ij (0) ~ sin0, similar to the planar tunnel 
junction. The effect of transparency for point contact between unconventional (d-wave) 
superconductors is studied in the papers [HI EHl EH1 133] ■ The non- locality of Josephson 
current in point contacts was investigated in |50] . The authors of [SU] concentrated on 
the influence of magnetic field on the zero voltage supercurrent through the junction. 
They found an periodic behavior in terms of magnetic flux and demonstrated that this 
anomalous behavior is a result of a non- locality supercurrent in the junction. This 
observation was explained theoretically in [51]. In this chapter we want to investigate 
theoretically the influence of electron reflection on dc Josephson effect in a ballistic 
point contact with transport current in the right and left banks which are separated 
by an interface (look at Fig0J). In Ref. [H2] for an ideal transparent point-contact in 
the impenetrable interface, it has been observed that at the phase differences close to 
the <p — 7r two antisymmetric vortex- like currents appear (see Fig|HJ). we want to study 
the effect of finite transparency (reflection) of the point-contact in the interface, on 
these vortex-like currents near the contact and at the phase difference <fi = n. We show 
that at low temperatures even a small reflection on the contact destroys the mentioned 
vortex-like current states, which can be restored by increasing of the temperature. In 
our system which will be investigated in this chapter, a point-contact between two 
massive superconductors (S-c-S junction) is considered as a possible Josephson weak 
links. For such systems it is convenient to use Kulik-Omelyanchouk method (SO] for 
the ballistic point-contact. The microscopic theory of the stationary Josephson effect 
in the ballistic point contacts between conventional superconductors was developed in 
[3*U] . Later, this theory was generalized for point-contacts between "d-wave" high-T c 
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supercondcutors in Ref. j!2]- The Josephson effect is a phase sensitive instrument for 
the analysis of an order parameter in novel (unconventional) superconductors, where 
current-phase dependencies Ij (0) may differ essentially from those in conventional 
superconductors [121 IE] • In some cases the model with ideal transparent point-contact 
does not correspond to the experiment totally, and the electron reflection should be 
taken into account. The influence of electron reflection on the Josephson current in 
ballistic point contacts was first considered by Zaitsev [31] . He had shown that reflection 
from the contact not only changes the critical value of current, but also the current- 
phase dependence Ij{4>) ~ sin (0/2) at low temperature which has been predicted 
in [30!- The current-phase dependence for small values of transparency, D <C 1, is 
transformed to the Ij (0) ~ sin 0, similar to the planar tunnel junction. In addition, 
the non-locality of Josephson current in point contacts was investigated in Ref. [50J. 
The authors of Ref. jSH| concentrated on the influence of magnetic field on the zero 
voltage supercurrent through the junction. They found an periodic behavior in terms 
of magnetic flux and demonstrated that this anomalous behavior is a result of a non- 
locality supercurrent in the junction. This observation was explained theoretically in 
[HUES EST Recently an influence of transport supercurrent, which flows in the contacted 
banks and is parallel to the interface, to the Josephson effect in point contacts has been 
analyzed theoretically[32 . It was found that a non-local mixing of two superconducting 
currents results in the appearance of two vortex-like current states in vicinity of the 
contact, when the external phase difference is ~ n. The Josephson current through 
superconducting weak link is a result of quantum interference between order parameters 
with phase difference 0. Obviously, the finite reflection R = 1 — D of electrons from 
the Josephson junction suppresses this interference and it must influence the vortex-like 
current states, which are predicted in [52*] . In this chapter, we study the effect of finite 
transparency on the current-phase dependence and distribution of the superconducting 
current near the ballistic point contact in the presence of homogeneous current states 
far from the contact. We show that at low temperatures (T — > 0) the electron reflection 
destroys the mentioned vortex-like current states even for a very small value of reflection 
coefficient flc 1. On the other hand we have found that, as the temperature increases 
the vortices are restored and they exist for transparency as low as D = | in the limit 
of T — > T c . The arrangement of the rest of this chapter is as follows. In Sec. fl3.2j) we 
describe the model of the point contact, quasiclassical equations for Green functions 
and boundary conditions. The analytical formulas for the Green functions are derived 
for a ballistic point contact with arbitrary transparency. In Sec. (j3.3|) we apply them 
to analyze a current state in the ballistic point contact. The influence of the transport 
current on the Josephson current and vice versa at the contact plane is considered. 
In Sec. ()3.4|) we present the numerical results for the distribution of the current in the 
vicinity of the contact. 
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3.2. Model and equations 

We consider the Josephson weak link as a microbridge between thin superconducting 
films of thickness d. The length L and width 2a of the microbridge, are assumed to be 
less than the coherence length £ - n the other hand, we assume that L and 2a are 
much larger than the Fermi wavelength and use the quasiclassical approach. There 
is a potential barrier in the contact, resulting in a finite probability for the electron that 
is to be reflected back. In the banks of superconductors a homogeneous current with a 
superconducting velocity v s flows parallel to the partition. We choose the z-axis along v s 
and the y-axis perpendicular to the boundary; y = is the boundary plane (see FigEJ). 
If the film thickness d <C £o then in the main approximation in terms of the parameter 
<i/£o the superconducting current depends on the coordinates in the plane of the film 
p = (y, z) only. The superconducting current in the quasiclassical approximation 

j(p, v s ) = AmeN(0)T £ (v f5i (v F) p, v s )) vp (53) 

is defined by the energy integrated Green matrix for the case of singlet superconductors 
which is following: 

w~ x ( 9l 92^2 \ /^x 

9{e, v f,P, v fl ) = .„ , (54) 

y io- 2 g 3 -9i J 

and for the case of the order parameter we have: 

A=f .° A ^ V (55) 
This Green matrix in the ballistic case satisfies the Eilenberger equations as follows 

[aU EH: 

f?—?^ + «'A*£ 2 (n) - iA n g 3{n) = 0; (56) 



dt 

~dt 



V^r- + 2e02(„) - 2iA n ^i (n) = 0; (57) 



- 2eg 3{n) + 2iA* n g l{n) = 0; (58) 

where, t = y/\v y \ on the Fermi surface, n = sgn(v y ) and n — 1, 2 label the left and right 
hand superconducting bulks, respectively. Using the quasiclassical approximation, we 
select the solution the for this problem as follows: 

9\(n) = £ + a n exp (-2sQt) ; (59) 
iA n 

92( n ) = + K exp (-2sQt) ; (60) 

93(n) = t -^ + d n exp(-2snt); (61) 

where, s = sgn(y) and Q = \Je 2 + |A| 2 . Here N(0) is the density of states at the Fermi 
level, e = e m + ipF-v s , vp and are the electron velocity and momentum on the Fermi 
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A exp(f) 



Figure 4. Model of the point-contact with finite reflection coefficient as a slit in the 
thin impenetrable insulating partition. 



surface, e m = (2m + 1)tvT are the Matsubara frequencies, m is an integer number, v s 
is the superfluid velocity and T is the temperature. Eqs. (|3.2|) should be supplemented 
by the equation for the superconducting order parameter A 

A(p, v„ T) = 2n\T £ (^(v F , p, v s )) Vf (62) 



m>0 



where A is the constant of pairing interaction and (••■) v is the averaging over directions 
of v F . After substitute in the Eilenberger equation ((TJ), we obtain: 



A n 
A* 

^ i 



A* 



exp (— 2sfl n t) ; 
exp (— 2sfi n t) ; 



(63) 
(64) 
(65) 



Q n \e + r]sQ n 

As it was shown in [SHj in the zero approximation in terms of the small parameter a/£o "C 
1 for a self-consistent solution of the problem it is not necessary to consider Eq. ([62)1 . 
The model, in which the order parameter is constant in the two half-spaces A(p, v s , T) = 
A(v s , T) exp(sgn(y) ^) ((f) is the phase difference between superconductors), can be used. 
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In the same approximation the velocity v s does not depend on the coordinates. 
The Eq. (|62)l enables us to calculate a spatial distribution of the order parameter A(p) 
in the next order approximation in terms of the parameter a/£o- Solutions of Eqs. (JT|) 
should satisfy Zaitsev's boundary conditions ([SI]) across the contact y — 0, \z\ < a 
and specular reflection condition for y — 0, \z\ > a. In addition, far from the contact, 
solutions should coincide with the bulk solutions. The Zaitsev boundary conditions 
[31] have been considered in [HI], but some improvements are necessary for using 
these boundary conditions. These improvements have been done in |55] . The Zaitsev 
boundary conditions at the contact can be written as [33 EH 155] 

d l = d r = d (66) 

= ds i2 (67) 



D 



D 



d^ 

:i+ 5 )« 
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where 

s r = g r (v F ,y = 0)+g r (v F ',y = 0) (68) 

d r = g r (v F ,y = 0)-g r (v F ',y = 0) (69) 

with vp' being the reflection of vp with respect to the boundary and D is the 
transparency coefficient of point contact. Indexes I and r denote that the Green function 
are taken at the left (y = — 0) or right (y = +0) hand from the barrier. Similar relations 
also hold for s' 1 and d . The first boundary condition implies that the antisymmetric 
part of Green function is continuous, this is a form of charge conservation, because 
the antisymmetric part of Green function is related to the current directly. But the 
second boundary condition means discontinuity in symmetric part of Green function. 
This discontinuity is the result of potential barrier as the interface. In general, D could 
be be momentum dependent. For simplicity in our calculations we assumed that D is 
independent of the Fermi velocity direction. 



3.3. Current-phase dependencies for Josephson and tangential currents 



Making use of the solution of Eilenberger equations (0), we obtain the following 
expression for the current density (J53j) at the slit [5^] : 



j(y = 0, \z\ <a)= 4neN(0)Tv F ^ / v F 



e£l — ir/DA 2 sin | cos 



m>0 



■~2 



+ A 2 



1 - £>(sin 




(70) 



where, Q 



-2 



A 2 , v = v F /v F is the unit vector and rj = sgn(v y ). In the case, 



v s ^ 0, the current (|7()jl has both jj and ] z components. The tangential current j z 
depends on the order parameters phase difference and is not equal to the transport 
current j F on the banks, in other words the total current is not equal to the vector sum 
of Josephson and transport currents. For the case v s = 0, at the contact the tangential 
current is zero and the normal component, i.e. the Josephson current is as found for the 
finite transparent contact in [33]. Detaching explicitly the Josephson current jj and the 
spatially homogeneous (transport) current jy that is produced by the superfluid velocity 
v a , we can write the current as the sum of three terms: jj, jr, and the "interference" 
current j» nt . Also we have 

j = jj(0,£>, V s 



+ Jt(v s ) +j int ((j),D,v s 



(71) 



The "interference" current takes place in the vicinity of the contact, where both coherent 
currents }j((f>) and jr(v s ) exist (see also the next subsection). At first we consider the 
current density ([70)) for temperatures close to the critical temperature (T c — T <^ T c ). 
From Eqs. (|70|) at the contact and for the temperatures close to the critical temperature 
we have: 



m - ir]DA 2 sin f cos f 



m>0 



e 2 



A 2 



1 - D(sin 



(e*) 2 + A 2 
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-£>(sinf 




(i*) 2 + A 2 


[i 


- £>(sinf 


) 2 ] 



(72) 
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Figure 6. Tangential current j z versus phase <j> for T/T c = 0.1 and q = 0.5. 

where, jo = 47r \e\ N(0)vpT c and numerator and denominator of the Green function 
fraction has been multiplied by the expression (e*) 2 + A 2 1 — D(sin |) , to escape of 
the complexity problems. It is well-known that for the temperatures close to the critical 
temperature, A(T — > T c ) « T c . Consequently, we use the Taylor expansion in terms 
of the small parameter A ^ T ~* Tc > . g a t the contact we have: 

-L c 

we obtain |56| : 



jr(v s ) = - 
j int ((f),D, v, 



-AD sin 0e y 



1 

-AA;e 2 
3 

1 



-AkD(l - cos<b)e z . 



(74) 
(75) 
(76) 
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Figure 7. Tangential current j z versus the transparency D at = ir and q = 0.5. 

where A = k = 14 j( 3 ) ej is the unit vector in the i— direction. This 

consideration shows how the current is affected by the interplay of Josephson and 
transport currents. At the contact the "interference" current j int is anti-parallel to 
jr and if the phase difference <ft = n, ji nt = —2D}t- When there is no phase difference 
(at (p = 0), we obtain j int = 0. So at transparency values D up to | the total tangential 
current at the contact flows in the opposite direction to the transport current. Thus, 
for such D in the vicinity of the contact, two vortices should exist. At arbitrary 
temperatures T < T c the current-phase relations can be analyzed numerically. In our 
calculations we define the parameter, q, in which q = and A = A(T = 0, v s — 0). 
The value of q can be in the range < q < q c and it's critical value q c , corresponds to 
the critical current in the homogeneous current state j2Zj- At T = 0, q c = 1 and the 
gap A does not depend on q. In FigEland Fig® we plot the Josephson and tangential 
currents at the contact as functions of at temperatures far from the critical (namely, 
T = 0.1T C ) and for q = 0.5 and for different values of transparency D. Far from <ft = n, 
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Figure 8. Vector plot of the current for <j) = w, q = 0.5, T/T c = 0.1 and D = 0.95. 
Axes are marked in units of the contact size a. 



the tangential current is not disturbed by the contact, it tends to its value on the 
bank. The Josephson current-phase relation is the same as when the transport current 
is absent. However, when tends to ir, for the highly transparent contact (D = 1,0.9) 
the tangential current becomes anti-parallel to the bulk current. But for D = 0.7 the 
"interference" current is strongly suppressed and the tangential current flows parallel to 
the bulk current. In FigEJ we plot j z (D) = jr+jint at = 7r for different temperatures. 
These plots show that by increasing the temperature a counter-flow j z (D) < exists in 
a wider interval of transparency D c (T) < D < 1 and D c (T — > T c ) — > |. This numerical 
result coincides with analytical results ()75I76|) . 
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Figure 9. Vector plot of the current for (j> = tt, q = 0.5, T/T c = 0.1 and D = 0.7. 



3.4- Spatial distribution of the current near the contact 

In this subsection we consider the spatial distribution of the current near the orifice. 
The superconducting current (JoHj) can be written as 



T 



j(p, v s ) = jo— v s )) 



V F 



(77) 



m>0 



where, j = 4ii\e\ N(0)vfT c . We should note that although the current (|77jl depends 
only on the coordinates in the film plane, the integration over velocity directions v is 
carried out over all of the Fermi sphere as in a bulk sample. This method of calculation is 
correct only for specular reflection from the film surfaces when there is no back scattering 
after electron interaction with them. At a point, p = (y,z), all ballistic trajectories 
can be categorized as transit and non-transit trajectories (see,Fig0J). For the transit 
trajectories "1" (their reflected counterparts marked by "3" in Fig0J) a projection vy of 
the vector v to the film plane belongs to the angle at which the slit is seen from the 
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Figure 10. Vector plot of the current for (j) = tt, q = 0.5, D = 0.7, and T/T c = 0.85. 

point p, vy G «(p), and for non-transit (marked by "2" in Fig0J Vn ^ a;(p)- For transit 
trajectories the Green functions satisfy boundary conditions on both banks and at the 
contact. The non-transit trajectories should satisfy the specular reflection condition [or 
Zaitsev's boundary conditions ()66 )) -([67 )) for D — at y — 0, \z\ > a]. Then for the 
current at T c — T <C T c we obtain an analytical formula [56J: 

j(p,(p,D,v 8 ) =j c D (sin 0vsgn(t> y ) + fc(l - cos0)vi£)~ -j c A; (vv y )~ (78) 

where, j c (T, v s ) = Elg!M2)EE a ^' v ^ . To illustrate how the current flows near the contact, 
we plot the FigJH] and FigEl for <fi = n and temperatures much smaller than critical 
(T/T c = 0.1), and for different values of transparency. At such value of the phase 
there is no Josephson current and at the large D = 0.95 the current is disturbed in such 
a way that there are two anti-symmetric vortices close to the orifice (see Fig. [HJ. For 
the such temperature at D = 0.7 the vortices are absent in Fig. El Near the critical 
temperature (T/T c = 0.85) the vortex-like currents are restored for D = 0.7 (see FigHUjl. 
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Figure 11. Tangential current j z versus the temperature T at <j> = it and q = 0.5. 



Far from the orifice (at the distances I ~ £o ^ a ) the Josephson current is spread out 
and the current is equal to its value at infinity. Considering the current distributions 
and current-phase diagrams, we observed that: 

1) . For fixed values of temperature and superfluid velocity, by decreasing the 
transparency the vortex-like current disappears at D = D c (T) ; 0.5 < D c (T) < 1 

2) . For intermediate values of transparency D (D C (T) < D < 1) by increasing the 
temperature the vortex-like currents, which were destroyed by the effect of electron 
reflection at the contact, may be restored. 

It is clear that both Josephson and "interference" currents are the result of the 
quantum interference between two coherent states. By decreasing the transparency 
the interference effect will be weaker and these two currents will decrease, while the 
transport current will remain constant. On the other hand, the presence of vortices 
depends on the result of competition between transport and "interference" current. 
Thus, by decreasing the transparency the tunneling and consequently the " interference" 
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current will decrease and vortices may be destroyed (|75I76|) . Similar to the case D = 1 
in at high values of transparency, the "interference" current can dominate the 
transport current and tangential current can be anti-parallel to the transport current, 
thus the vortices appear. But for low transparency the tangential current will be parallel 
to the transport current and the vortices disappear. 

The second point is an anomalous temperature behavior of the effect. The vortices 
are the result of the coherent current mixing. One could expect that by increasing 
the temperature the vortices would disappear whereas, for intermediate values of 
transparency, by increasing the temperature the vortices will be restored. As considered 
in FigEland FigElfor the transparency D = 0.7 the vortices at low temperature are 
absent but at high temperature they are present. In the plots for tangential current 
versus transparency, Fig. 0we can observe this phenomenon (appearance of the counter- 
flow near the contact at high temperatures). 

Usually superconducting currents are monotonic and descendant functions of 
temperature. Josephson and transport currents have this property, but about the 
tangential current j y , the situation is totally different. At high values of transparency 
the j y has similar behavior to the two other currents, but at low and intermediate 
values of transparency at = 7r it has a non-monotonic dependence on the temperature 
and this is the origin of the anomalous temperature behavior of vortices. As the 
temperature increases, the tangential current first increases and then decreases. In 
FigE]we plotted the tangential current ( "interference"-!- transport current) versus the 
temperature for different values of transparency. We observed that for intermediate 
values of transparency 0.5 < D < 1, at low temperatures and = tt the tangential 
current has anomalous dependence on the temperature. The reason for this dependence 
is that the "interference" current flows in the opposite direction to the transport current. 
This current is suppressed by the reflection, but with increasing of the temperature it 
decreases slowly than the transport current. As a consequence of that with increasing 
of T the tangential current can change its sign and vortices appear. We found that for 
low values of transparency < D < 0.5 the "interference" current cannot dominate the 
transport current and in addition the tangential current has the same direction as the 
transport current for any temperature T < T c . 
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4. Weak link between unitary /-wave superconductors 

4-1. Introduction 

In this chapter the spin current in the Josephson junction as a weak-link (interface) 
between misorientated triplet superconductors will be investigated theoretically for the 
models of the order parameter in UPt^. Green functions of the system will be obtained 
from the quasiclassical Eilenberger equations. The analytical results for the charge and 
spin currents will be illustrated by numerical calculations for the certain misorientation 
angles of gap vector of superconductors. 

Triplet superconductivity has become one of the most interesting topics of condensed 
matter physics [SJ [7j, particularly in view of the recently discovered ferromagnetic 
superconductivity ^31 C3! The mechanism of pairing, physics of interaction and gap 
structure in this type of superconductors have been the subject of many experimental 
and theoretical works [5B1EH]. The Cooper pairing in the triplet superconductors has 
been investigated, for example, using the thermal conductivity in papers [HOI EE] and 
Knight shift experiments in papers jH2 Also, the Josephson effect in the point 
contact between triplet superconductors has been studied in paper These weak-link 
structures have been used to demonstrate the order parameter symmetry in Ref. [64J. 
Eventually, the /-wave symmetry of order parameter has been proposed for UPt% and 
SrzRuOi compounds. In addition, the spin polarized transport through the systems 
consisting of superconductors, normal metals, ferromagnetic layers and other structures 
as one of the modern topics of mesoscopic physics, has attracted much attention recently 
[Hill EH1 E3 EE]- In this chapter, the ballistic Josephson weak-link as the interface between 
two bulk of /-wave superconductors with different orientations of the crystallographic 
axes has been investigated. It is shown that the current-phase dependencies are totally 
different from the current-phase dependencies of the junction between conventional (s- 
wave) superconductors |30j and high T c (ci-wave) superconductors [IB] . It is found that 
for the certain values of the misorientation, the spin-current in the both directions, 
tangential and perpendicular to the interface, may exist and it has totally unusual 
dependence on the external phase difference. The effect of misorientation on the spin 
current is investigated. It is observed that the misorientation between gap vectors is the 
origin of the spin current. As the important result of this chapter, it is obtained that, 
at some of certain values of phase difference, at which the charge current is zero, the 
spin current has the finite value. Another result of this chapter is the capability of this 
proposed experiment for polarization of the spin transport using the junction between 
/-wave superconductors. Eventually, one of the states and geometries of our system can 
be used as a switch which is able to divide the spin and charge currents into two parts: 
parallel and perpendicular to the interface. 

The arrangement of the rest of this chapter is as follows. In Sec. (j4.2j) we describe 
our configuration, which is investigated. For a non-self-consistent model of the order 
parameter, the quasiclassiacl Eilenberger equations [HI] are solved and suitable Green 
functions are obtained analytically. In Sec. ([4.3)1 the obtained formulas for the Green 



Charge and Spin Transport in Superconducting Weak Links 



34 




Figure 12. Scheme of a flat interface between two superconducting bulks, which are 
misorientated as much as a. 



functions are used for calculation the charge and spin current densities at the interface. 
An analysis of numerical results will be done in Sec. ([4.4)1 . 

4-2. Basic Equations 

We consider a model of a flat interface y = between two misorientated /—wave 
superconducting half-spaces (FigH2J) as a ballistic Josephson junction. In the 
quasiclassical ballistic approach, in order to calculate the charge and spin current, we use 
"transport-like" equations [3T] for the energy integrated Green functions g(vp,r, e m ) 



v F V<? + 



£ m a 3 + iA,g 



0, (79) 
and the normalization condition 

99 = 1, (80) 

where e m = 7iT(2m + 1) are discrete Matsubara energies m = 0,1,2..., T is the 
temperature and is the Fermi velocity and 03 = (T3 <g> / in which &j (j = 1, 2, 3) 
are Pauli matrices. The Matsubara propagator g can be written in the standard form: 



9 



gi + gi<7 {g 2 + g 2 <x) i&2 
i&2 (93 + g3<5") g^ - d- 2 g40-a 2 
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where, the matrix structure of the off-diagonal self energy A in the Nambu space is 

A=f ° . d 7 2 V (82) 

In this chapter, the unitary states, for which dx d* = 0, is investigated. Also, the 
unitary states vectors di^ can be written as 

d„ = A n exp^ n , (83) 

where A 12 are the real vectors in the left and right sides of the junction. The gap (order 
parameter) vector d has to be determined from the self-consistency equation, near the 
Fermi surface: 

d (v F , r) = ttTN (0) £ (V (v F , v' F ) g 2 (v' F , r, e m )) (84) 

m 

where V (v>, v F ), is a potential of pairing interaction, (...) stands for averaging over the 
directions of an electron momentum on the Fermi surface and N (0) is the electron 
density of states at the Fermi level of energy. Solutions to Eqs. (fTHJ) and (jHljl must satisfy 
the conditions for Green functions and vector d in the bulks of the superconductors far 
from the interface as follow: 

g(±oo) = =^=- (85) 

/ ^+|d 2 ,i| 2 



d(±oo) =d 2il (v F )exp( T ^), (86) 



2 

where <fr is the external phase difference between the order parameters of the bulks. Eqs. 
(J79J) and (JS4"j) have to be supplemented by the continuity conditions at the interface 
between superconductors. For all quasiparticle trajectories, the Green functions satisfy 
the boundary conditions both in the right and left bulks as well as at the interface. 
The set of equations (|79J) and (|H4*|) can be solved only numerically. For unconventional 
superconductors such solution requires the information of the function ^(v F ,Vp). 
This information, as that of the nature of unconventional superconductivity in novel 
compounds, in most cases is unknown. Usually, the spatial variation of the order 
parameter and its dependence on the momentum direction can be separated in the form 
of A(vp, y) = A(y F )^/(y). It has been shown that the absolute value of a self-consistent 
order parameter and ^(y) are suppressed near the interface and at the distances of the 
order of the coherence length, while its dependence on the direction in the momentum 
space (A(vp)) remains unaltered {Hj. Consequently, this suppression doesn't influence 
the Josephson effect drastically. This suppression of the order parameter keeps the 
current-phase dependence unchanged but, it changes the amplitude value of the current. 
For example, it has been verified in Ref. [IB] for the junction between unconventional 
(i-wave, in Ref. jH] for the case of "/-wave" superconductors and in Refs. [60J [70] 
for pinholes in 3 He that, there is a good qualitative agreement between self-consistent 
and non-self-consistent results. Also, it has been observed that the results of the non- 
self-consistent investigation of ferromagnet-d-wave proximity structure in Ref. [ZTj are 
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Figure 13. Charge and spin current (s z ) versus the phase difference <j> for the planar 
state (|1U3II . geometry (i) and the different misorientations. Currents are given in units 
of j 0>e = |eiV(0)w F A (0) and j , s = jhN(0)v F A o (0) respectively. 



coincident with the experimental results of the paper [72] and the results of the non-self- 
consistent model in paper [33] are similar to the superffuid weak-link experiment [73]- 
In the paper [7T], they have investigated the proximity effect between a ferromagnet 
and a high— T c superconductor. They have solved the Eilenberger equation and using 
the obtained Green function investigated the Andreev bound states. The density of 
states in this system has been studied and the spatial oscillations have been observed 
in this nonselfconsistent paper. The results of this paper had been observed before in 
an experimental report in paper |72j . In addition, there are many published papers 
[HH1 EH [7T] and [7U [731 113 EI] in which, such approximation has been used for different 
systems containing unconventional superconductors and important analytical results 
have been obtained. In Refs. [^3 EH HI] Eilenberger equation has been solved and 
Bogoliobov- Degennes equation has been considered in papers [ZU UJ3 EH E7J non- 
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Figure 14. Charge and spin current (s y ) versus the phase difference <j> for the axial 
state (|102|) . geometry (ii) and the different misorientations (y-component). 



self-consistently. Consequently, despite the fact that this estimation cannot be applied 
directly for a quantitative analyze of the real experiment, only a qualitative comparison 
of calculated and experimental current-phase relations is possible. In our calculations, 
a simple model of the constant order parameter up to the interface is considered and 
the pair breaking and the scattering on the interface are ignored. We believe that under 
these strong assumptions our results describe the real situation qualitatively. In the 
framework of such model, the analytical expressions for the charge and spin current 
can be obtained for an arbitrary form of the order parameter. Also, we have done 
our calculations for the small misorientations a = ^j, a = ^ and a = jL. As the 
results of paper |2H], for the small misorientations, selfconsistent and nonselfconsistent 
calculations have the same results approximately. Consequently, authors of paper jlH] 
concluded that the nonselfconsistent formalism can be used for the junction between 
unconventional superconducting bulks with small misorientations. 
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4-3. Analytical Green functions 

The solution of Eqs. (fT9|) and (J84|) allows us to calculate the charge and spin current 
densities. The expression for the charge current is: 

j e (r) = 2meTN (0) ^ (v f9i (v F , r, e m )) (87) 

and for the spin current we have: 

j Si (r) = 2m{^)TN (0) ^ (v F (e igl (v>, r, e m ))) (88) 

m 

where, ej=(x, y,z). We assume that the order parameter does not depend on 
coordinates and in each half-space it equals to its value (jll2j) far from the interface in 
the left or right bulks. For such a model, the current-phase dependence of a Josephson 
junction can be calculated analytically. It enables us to analyze the main features of 
current-phase dependence for the different models of the order parameter of "/-wave" 
superconductivity. The Eilenberger equations (|79|) for Green functions g, which are 
supplemented by the condition of continuity of solutions across the interface, at y — 0, 
and the boundary conditions at the bulks, should be solved for a non-self-consistent 
model of the order parameter analytically. In a ballistic case the system of equations 
for functions gi and gj can be decomposed on independent blocks of equations. The set 
of equations which enables us to find the Green function gi is: 

v F kVgi = i (d ■ g 3 - d* • g 2 ) ; (89) 

v F kS7g- = -2(dx fl3 + d*xg 2 ) ; (90) 

v F kVg 2 = -2e m g 2 + 2igid + d x g_; (91) 

v F kVg 3 = 2e m g 3 - 2i 9l d* + d* x g_; (92) 

where g_ = gi — g 4 . The Eqs. (}9"2j) can be solved by integrating over ballistic 

trajectories of electrons in the right and left half-spaces. The general solution satisfying 
the boundary conditions (jl()9j) at infinity is 

^(n) _ £m + ^ (_2sQ n t) ; (93) 
g W = C n exp(-2^ n t); (94) 

g f } = ^ - 2m : d " + _ W 9 XCn ex p (- 2sn ^ ; ( 95 ) 

where t is time of flight along the trajectory, sgn (t) = sgn (y) = s and r] = sgn (v y ) . 
By matching the solutions (j93H96j) at the interface (y — 0, t — 0), we find constants a n 
and C n . Indices n = 1, 2 label the left and right half-spaces respectively. The function 
9i (0) = ( — 0) = g^ (+0) , which is a diagonal term of Green matrix and determines 
the current density at the interface, y = 0, is as follows: Two diagonal terms of Green 
matrix which determine the current densities at the interface, y = 0, are following. For 
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Figure 15. Charge and spin current (s y ) versus the phase difference <j> for the axial 
state l|lU2f) . geometry (ii) and the different misorientations (x-component) 



the relative term to the charge current we obtain: 

( , _ gmffli + ^2) cos/3 + vq sin j3^h^h + g|J , ^ 

^ U ~ 177 sin /fc^ + n 2 ) + cos/?(fi 1 fi 2 + e 2 J + A,A 2 1 j 

and for the case of spin current we have: 

gi (0) = (98) 
M[(B - l) 2 exp(^)(rA + e m )(v^2 + e m ) - (B + l) 2 exp(-z/3)(^ 2 - £ m )(rA - e m )] 



where 77 = s#n (u w ), fi n = \J e 2 m + |d n | , f3 = ^ - ^2 + 0, 

^ _ r?£ m (Qi + fl 2 ) cos/3 + isin/3(^if2 2 + e 2 m ) 
~ ir) sin pemitti + fi 2 ) + cos/3(fiifi 2 + e*J + A X A 2 ' 

^ = A^QB-ljexp^) + A 1 A 2 (B + 1) expj-iP) 
(rjQi - e m ){i]Q 2 - e m ) (77^1 + e m ) (r?fi 2 + e m ) 
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Figure 16. Charge and spin current (s y ) versus the phase difference <j> for the axial 
state H1U2|) . geometry (ii) and the different misorientations (^-component) 



and 

M = '— 1 = 101 

(A + 2 J B)|d 1 | 2 |d 2 | 2 

Also, n — 1, 2 label the left and right half-spaces respectively. We consider a rotation 
R only in the right superconductor (see, FigH2J), (i.e., d 2 (k) = -Rdi(_R _1 k), k is the 
unit vector in the momentum space). The crystallographic c-axis in the left half-space 
is selected parallel to the partition between the superconductors (along z-axis in FigEJ). 
To illustrate the results obtained by computing the formula (I97I98|) . we plot the current- 
phase diagrams for the different models of the "/-wave" pairing symmetry (jl02H03|) and 
for two different geometries. These two geometries are corresponding to the different 
orientations of the crystals in the right and left sides of the interface (see, Fig IT2"|) : 

(i) The basal afo-plane in the right side has been rotated around the c-axis by a; Ci||c2. 

(ii) The c-axis in the right side has been rotated around the axis perpendicular to the 
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interface (y-axis in Fig|T^|) by a; bi||b 2 . 

Further calculations require a certain model of the gap vector (vector of order parameter) 
d. 

4-4- Numerical results 

In this chapter, two most possible forms of the /-wave order parameter vector in UPts 
are considered. The first model which is successful to explain the properties of the B- 
phase of UPt s is the axial state. This sate describes the strong spin-orbital coupling 
with vector d directed along the c axis of the lattice and it is: 

d(T, vj,) = A (T)zk z (k x + iky) 2 . (102) 

The coordinate axes x, y, z here and below are chosen along the crystallographic axes 
a, b, c in the left side of FigHJl The function A = A (T) in Eq. {TPS} and below 
describes the dependence of the order parameter d on the temperature T (our numerical 
calculations have been done at the temperatures close to the T = 0). The second model 
of the order parameter which describes the weak spin-orbital coupling in UPt^ states, 
is the unitary planar state. The planar model of gap vector is: 

d(T, v F ) = A (T)k z (± (k 2 x - k 2 y ) + y2k x k y ). (103) 

Using these two models of order parameters (jl02H03|) and solutions to the Eilenberger 
equations (JiTTj) and (J98|) . we have calculated the spin current and charge current densities 
at the interface numerically. These numerical results are listed below: 

1) The spin current can be present, only when misorientation between gap vectors 
exists. Because in our Green function the spin current is proportional to the 
"cross product" between the left and right gap vectors. For instance, the spin current 
for the case of the axial state (J102j) and geometry (i) is zero, because both of the gap 
vectors are in the same direction (z). (Geometry (i) is a rotation as much as a, around 
the z axis). 

2) In FigEHit is shown that for the planar state and geometry (i), it is possible to observe 
the current of s z in the direction perpendicular to the interface, but in Figs J14lIo1 and 
HEl it is demonstrated that, for the axial state and geometry (ii), only the current of s y 
can be observed. Consequently, this kind of junction can be applied as a polarizer or 
filter for the spin currents. 

However, for the planar state and geometry (ii), all terms of the spin current (s x , s y and 
s z ) can be observed (see Eql^Hj). 

3) In Figs ll4ll2"Ul (planar states), it is shown that the value of the phase differences in 
which the currents are in the maxima, minima and zero values, are not very sensitive 
to the misorientation angle a, while the amplitude of maxima and minima, are strongly 
dependent on the value of misorientation a. 

4) In the Figs ll3pi4l and FigHT} ITS} while the charge currents are the odd functions 
of with respect to the line of = 7r, the spin currents are even functions of 
the phase difference; j e (0 = tt + 5<p) = — j e (0 = 7r — 5(f)) and for the spin current 
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Figure 17. Tangential spin current (s z ) versus the phase difference 4> for the 
planar state (|l(J3f) . geometry (ii) and the different misorientations. The perpendicular 
component (y-direction) of the spin current is absent. 



= 7r + 5<p) = j Sl (0 = 7T — S(f>). On the contrary, in the FigslTBI and UTA the charge 
and spin currents are even and odd functions of with respect to the line of = tt, 
respectively; j e (0 = tt + 5<p) = j e (0 = tt - 5<p) and j Si (4> = tt + 5<p) = = 7r - 50). 

5) In Fig JTHl the perpendicular component of the spin and charge current in terms of 
the external phase difference for the case of the planar state (jlOHjl . geometry (i) and 
for two different misorientations are plotted. The solid line is the charge current-phase 
dependence Also, at the = 0, = tt and = 2tt, the charge current (Josephson 
current) is zero while the spin current has the finite value. 

6) The perpendicular component of the charge (Josephson current) and spin current 
for the case of the axial state (jl()2j) and geometry (ii) are plotted in Figll4[ and the 
tangential components of them, are plotted in Figs J15li7fl The charge current-phase 
diagrams have been obtained before in paper jH] and they are totally different from the 
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Figure 18. The perpendicular component of the charge current (Joscphson) versus 
the phase difference <j> for the planar state (|il)3[) . geometry (ii) and the different 
misorientations. The tangential components (x and z-directions) are absent. 

case of conventional superconductors in the paper [30]. At the phase values of <fi = 0, 
= 7r and = 2tt, in which the charge current is exactly zero, the spin current has the 
finite values and may select its maximum value. In Figs ll3JiH and specially FigEOl for 
a small value of misorientation we have a very long but narrow peak in the spin current 
phase diagram, close to the = n. 

7) Both the planar state with geometry (i) and the axial state with geometry (ii) can 
be applied as the filter for polarization of the spin transport (see Fig dlHllTKJ) . the former 
transports only the s z but the latter case flows the s y (see 1E)S|) . In addition, the planar 
states with geometry (ii) can be used as a switch for the spin and charge current into 
two directions: parallel and perpendicular to the interface. In this case, the spin and 
charge currents select only one of the directions parallel or perpendicular to the interface. 
Namely, it is impossible to observe the tangential and perpendicular components of the 
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Figure 19. Tangential spin currents (s x ) versus the phase difference 4> for the 
planar state (|l(J3f) . geometry (ii) and the different misorientations. The perpendicular 
component (y-direction) of the spin current is absent. 

currents at the same time for planar state with geometry (ii)(FigsEHHI|)- 
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Figure 20. Perpendicular component of the spin current (s y ) versus the phase 
difference <f> for the planar state (|ll)3fl . geometry (ii) and the different misorientations. 
The tangential components (x and z directions) are absent. 



5. Weak link between PrOs^Sb^ supercondcuting banks 

5.1. Introduction 

The "(p + h)— wave" form of pairing symmetry has been considered for the 
superconductivity in PrOs^Sb^ compound, recently. In this chapter, a stationary 
Josephson junction as a weak-link between PrOs^Sb^ triplet superconductors is 
theoretically investigated. The Eilenberger equation is solved for two distinct models of 
the order parameter (A and B— phases). The spin and charge current-phase diagrams are 
plotted and the effect of misorientation between crystals on the Josephson, spontaneous 
and spin current is studied. It is obtained that such experimental investigations of the 
current-phase diagrams can be used to test the pairing symmetry in the above-mentioned 
superconductors. In this chapter, it is shown that this apparatus can be applied as a 
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polarizer for the spin current. Furthermore, it is observed that at some certain values 
of the phase differences in which the charge current is zero, the spin current exists while 
carriers of both of charge and spin are the same (electrons) . 

The pairing symmetry of the recently discovered superconductor compound PrOs^Sbu 
is an interesting topic of research in the field of superconductivity (HI ESI EBj- 
Superconductivity in this compound has been discovered in papers [7H1 [HOI EJ] and 
two different phases (A and B) have been considered for this kind of superconductor in 
papers [251 EH]- Although people at first have considered the spin-singlet u (s+g)— wave" 
pairing symmetry to this superconductor [25J but later it has been specified that the spin- 
triplet is the real pairing symmetry of PrOs^Sbu complex [SJIH2!- Authors of paper [82J, 
using the knight shift in NMR measurement estimated the spin-triplet pairing symmetry 
for the superconductivity in PrOs^Sb^. Consequently, the "(p + h)— wave" proposed 
for the pairing symmetry of the superconductivity in PrOs4,Sbi2 compound, recently 
jH]. In this chapter, the self-consistent equation for the gap vectors (BCS gap equation) 
have been solved for the finite temperature numerically and for the temperatures close 
to the zero and the critical temperature analytically. For this compound and using the 
u (p+h)— wave" symmetry for the order parameter vector (gap function), the value of the 
order parameter at the zero temperature has been obtained for both A and -B-phases, 
in terms of the critical temperature of this type of superconductivity (T c ). In addition, 
the gap dependence on the temperature close to the zero and critical temperatures have 
been obtained. Authors of paper jH] have investigated the critical magnetic field and 
the temperature dependence of critical field, specific heat and heat conduction. 
Also, the Josephson effect in the point contact between triplet superconductors has 
been studied in paper p. In this chapter, the effect of misorientation on the charge 
transport has been studied and a spontaneous current tangential to the interface between 
the /—wave superconductors has been observed. Additionally, the spin-current in the 
weak-link between the /—wave superconductors has been investigated in the paper 
|H3j . The authors of paper jHH], have proposed this kind of weak-link device as the 
filter for polarization of the spin-current. These weak-link structures have been used to 
demonstrate the order parameter symmetry in paper 

In this chapter, the ballistic Josephson weak-link via an interface between two bulk 
of "(p + h)— wave" superconductors with different orientations of the crystallographic 
axes is investigated. It is shown that the spin and charge current-phase diagrams are 
totally different from the current-phase diagrams of the junction between conventional 
(s-wave) superconductors [30J, high T c (d-wave) superconductors [48J and from the spin- 
current phase diagrams in the weak-link between the /-wave superconductors [83J. In 
this weak-link structure between the "(p + h)— wave" superconductors, the spontaneous 
current parallel to the interface as the characteristic of unconventional superconductivity 
can be present. The effect of misorientation on the spontaneous, Josephson and spin 
currents for the different models of the paring symmetry (A— and B— phases in Fig J5.1|) 
are investigated. It is possible to find the value of the phase difference in which the 
Josephson current is zero but the spontaneous current tangential to the interface, which 
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Figure 21. A— phase (left-above), B— phase (right-above) order parameters, direction 
of a, b, c and z unit vectors (left-below) and chemical structure of PrOsiSb 12 
molecule [25]. A— and B— phases are high field (high temperature) and low field 
(low temperature) phases, respectively [2*7) . 

is produced by the interface, exists. In some of configurations and at the zero phase 
difference, the Josephson current is not zero and it has a finite value, this finite value 
corresponds to a spontaneous phase difference which is related to the misorientation 
between the gap vectors. 

Finally, It is observed that at the some certain values of the phase differences <fi in which 
the charge current is zero the spin current exists and vise versa. In addition, in this 
configuration in which the gap vectors are selected along the z direction and the unit 
vector perpendicular to the interface is y direction, only the spin current of the s y can 
be present and the other terms of the spin current are absent totally. Consequently, this 
structure can be used as a filter for polarization of the spin transport. Furthermore, our 
analytical and numerical calculations have shown that the misorientation is the origin 
of the spin current and in the absence of the misorientation spin current is absent while 
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Figure 22. Perpendicular component of current (Josephson current) versus the 
phase difference <j> for A and B— phases, geometry (i), T/T c = 0.08 and the different 
misorientations. Currents are given in units of jo = ^eN(0)vp^a(0)- 



the charge current flows. 

In Sec. (j5.2j) the obtained formulas for the Green functions will be discussed and an 
analysis of numerical results will be done. 

5.2. Discussions 

In this chapter again we consider a model of a flat interface y = between two 
misorientated "(p + h)— wave" superconducting half-spaces fFiglT^j) as a ballistic 
Josephson junction. As same as the previous chapter, in our calculations, a simple 
model of the constant order parameter up to the interface is considered and the pair 
breaking and the scattering on the interface are ignored. The solution of Eqs. (JIJ and 
()84|) allows us to calculate the charge and spin current densities. Again, we assume that 
the order parameter does not depend on coordinates and in each half-space it equals 
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Figure 23. Perpendicular component of current (Josephson current) versus the 
phase difference <f> for A and B— phases, geometry (ii), T/T c = 0.08 and the different 
misorientations . 

to its value (|112j) far from the interface in the left or right bulks. For such a model, 
the current-phase dependence of a Josephson junction can be calculated analytically. It 
enables us to analyze the main features of current-phase dependence for the different 
models of the order parameter of "(p + h)— wave" superconductivity To illustrate the 
results obtained by computing the formula (|97j). we plot the current-phase diagrams 
for the different models of the "(p+ h)— wave" pairing symmetry (JKJ4I1U5)) and for two 
different geometries. These geometries are corresponding to the different orientations of 
the crystals in the right and left sides of the interface (FigJT2|): 

(i) The basal afe-plane in the right side has been rotated around the c-axis by a; Cx || c 2 - 

(ii) The c-axis in the right side has been rotated around the 6-axis by a (y-axis in 
FigED; bi||b 2 . 

Further calculations require a certain model of the gap vector (vector of order parameter) 
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Figure 24. Tangential (to the interface) component of current versus the phase 
difference <fi for A and B— phases, geometry (ii), T/T c = 0.08 and the different 
misoricntations(a;— component). 



d. 

In this chapter, two forms of the "(p + h)— wave" unitary gap vector in PrOs^Sb^ 
are considered. The first model to explain the properties of the A-phase of PrOs^Sbw 
is: 

d(T, k) = A (T)(k x + ik y )\{\ -kt-kt- kt)z (104) 

The coordinate axes x, y, z here and below are chosen along the crystallographic axes 
a, b, c in the left side of FigUJ The function A = A (T) in Eq. (ITM|l and below 
describes the dependence of the gap vector d on the temperature T. The second model 
to describe the gap vector of the S-phase of PrOs^Sbu is: 

d(T,k) = A (T)(k x + ik y )(l-k A z )z (105) 

Our numerical calculations are done at the low temperatures, T/T c = 0.08, and 
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Figure 25. Tangential (to the interface) component of current versus the phase 
difference <fi for the both A and B— phases, geometry (ii), T/T c = 0.08 and the different 
misorient ations (z — component ) . 



we have used the formulas ln(A(T)/A(0)) = — (^y) 3 for the A— phase and 

ln(A(T)/A(0)) = -^f^^o)) 3 for the 5-phase, from the paper j8], for temperature 
dependence of the gap functions Ao(T) at the low temperatures (T << T c ). Also, in the 
mentioned paper [H] the value of A (T) in terms of the critical temperatures has been 
calculated for both A and B— phases. They are A(0)/T c = 2.34 and A(0)/T c = 1.93 for 
A and B— phases, respectively. Using these two models of order parameters f|lU4llU5|) 
and solution to the Eilenberger equations ()97|98|) . we have calculated the current density 
at the interface numerically. These numerical results are listed below [84J: 
1) In Fig |22l the perpendicular component of current (perpendicular to the interface) 
which is called Josephson current is plotted for both A and B— phases, geometry (i), 
misorientations a = it/ '4 and a = ir/6. It is observed that the critical values of 
current for B— phase is greater than A— phase. Also, in spite of junction between the 
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Figure 26. Perpendicular component of current (Josephson current) versus the phase 
difference <p f° r geometry (i), jr — 0.08, a = j, A— phase, B— phase and p— wave 
pairing symmetry. 



conventional superconductors and planar Josephson junction, at the = 0, the current 
is not zero. The current is zero at the phase difference value = 0o, which depends 
on the misorientation between the gap vectors. In Fig|221 the value of the spontaneous 
phase difference 0o is close to the misorientation a. 

2) In Fig 1221 Josephson current is plotted for both A and B— phases, geometry (ii) and 
different misorientations. Again, the maximum value of the current for the B— phase 
is greater than A— phase. Increasing the misorientation between the gap vectors, the 
maximum value of current decreases. It is demonstrated that at the phase difference 
values = 0, = 7r and = 2tt, the Josephson current is zero while, both of spontaneous 
and spin currents are not zero and have the finite value. Increasing the misorientation 
between the gap vectors decreases the derivative of the current with respect to the phase 
difference (4r) close to the « 



7T. 
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Figure 27. Tangential component of spin (s y ) current (x— component) versus the 
phase difference <f> for geometry (ii), ^- — 0.08 and different misorientations between 
the A and B— phase of "(p + h)— wave" pairing symmetry. 



3) In Figs l24l2B*l the tangential components of charge current (x and z— components) 
in terms of the phase difference are plotted. It is seen that at the = 0, = ir and 
= 2tt in which the Josephson current is zero the parallel spontaneous currents have 
the finite values. Although the perpendicular component of charge current (see FigEBJ) 
is an odd function of the phase difference with respect to the line of = it but the 
parallel charge currents are even functions at the phase differences close to = n (see 
FigEland FigESJ- 

4) In Fig|2nj the Josephson current in terms of the phase difference is plotted for the case 
of p— wave, and A and B— phases of u (p + h)— wave". The p— wave pairing symmetry 
as the first candidate for the superconducting state in Sr 2 Ru0 4i is as follows [T§] : 

d(T,k) = A Q (T)(k x + ik y )z (106) 

It is observed that the maximum value of Josephson current (j y ) of junction between 
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Figure 28. Tangential component of spin (s y ) current (z— component) versus the 
phase difference <f> for geometry (ii), ^- — 0.08 and different misorientations between 
the A and B— phase of "(p + h)— wave" pairing symmetry. 



the p— wave superconductors, is greater than the B— phase of "(p + h)— wave" and the 
Josephson current of second is greater than its A— phase counterpart. Also, the place 
of the zero of the current is at the spontaneous phase difference which is close to the 
misorientation <fto = a (look at the FigJ32|l. 

5) In FigsEZEHl the tangential spin (s y ) currents in terms of the phase difference is 
plotted for different misorientations and geometry (ii). By increasing the misorientation 
the maximum value of spin current increases. In spite of the charge current for this 
state, the spin current at the phase differences (f> = 0, (f> = tt and (p = 2tt is exactly zero 
(compare Figs|21 and ESI with Figsl27l and ¥M . 

6) In FigsEHEOl the perpendicular component of charge and spin current (j y and j sy ) are 
plotted for different misorientations and A and B— phases respectively. An interesting 
case in our observations is the finite value of the perpendicular spin current at the = 0, 
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Figure 29. Perpendicular component of charge and spin current (j y and j sy ) versus 
the phase difference <j) for A— phase, geometry (ii), ^- = 0.08, a = ? and a = j. 

<p = it and = 2tt at which the perpendicular charge current (j y ) is zero (see Figs 1^1 
andlSOl). 

7) In FiglHTl the perpendicular component of the spin current is plotted for p— wave, 
A and B— phases of u (p + h)— wave" pairing symmetries and for a specified value of 
misorientation a = In both FigsEHlandl^the maximum value of current of junction 
between the p— wave is greater than B— phase and B— phase has the maximum value 
greater than junction between the A— phase. This different character of the current- 
phase diagrams enables us to distinguish between the three states. Also, It is observed 
that at the phase differences <fi = 0, = 7r and = 2tt, the spin current has the finite 
value and may have its maximum value. This is a counterpart of Fig|2Hl in that figure 
the charge currents are zero at the mentioned value of the phase difference but the spin 
current has the finite value. 

Furthermore, our analytical and numerical calculations have shown that the origin of 
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Figure 30. Perpendicular component of charge and spin current (j y and j sy ) versus 
the phase difference 4> f° r B— phase, geometry (ii), jr — 0.08, a = ^ and a = j. 

the spin current is misorientation between the gap vectors (cross product in Eq 198(1. 
Thus the spin current in the weak link between geometry (i) misorientated crystals 
is zero. Because the geometry (i) is a rotation by a around the z— axis and both of 
the left and right gap vectors are in the same direction and cross product between 
them is zero. Also, it is shown that in this structure (y— direction is perpendicular to 
the interface, gap vectors c— axis are selected in the z— direction and rotation is done 
around the y— direction) only the current of the s y flows and other terms of the spin 
current are totally absent. So, this kind of weak-link experiment can be used as the 
filter for polarization of spin transport. Since the spin is a vector, the spin current is a 
tensor and we have the current of spin s y in the three x, y and z directions. 
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Figure 31. Perpendicular component of spin current (j ay ) versus the phase difference 
for geometry (ii), ^- — 0.08, a = j, A— phase, B— phase and p— wave pairing 
symmetry. 

6. Weak link between nonunitary triplet superconductors 

6.1. Introduction 

A stationary Josephson effect in a weak-link between misorientated nonunitary triplet 
superconductors is investigated theoretically. The non-self-consistent quasiclassical 
Eilenberger equation for this system has been solved analytically. As an application 
of this analytical calculation, the current-phase diagrams are plotted for the junction 
between two nonunitary bipolar /—wave superconducting banks. A spontaneous current 
parallel to the interface between superconductors has been observed. Also, the effect of 
misorientation between crystals on the Josephson and spontaneous currents is studied. 
Such experimental investigations of the current-phase diagrams can be used to test the 
pairing symmetry in the above-mentioned superconductors. 
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In the recent years, the triplet superconductivity has become one of the modern 
subjects for researchers in the field of superconductivity [2J 13 EH]- Particularly, the 
nonunitary spin triplet state in which, Cooper pairs may carry a finite averaged intrinsic 
spin momentum, has attracted much attention in the last decade [T51 EE] • A triplet 
state in the momentum inverse space (k) can be described by the order parameter 
A(k) = i(d(k) ■ o)(j y in a 2x2 matrix form in which &jS are 2x2 Pauli matrices. The 
three dimensional complex vector d(k) demonstrates the triplet pairing state exactly 
In the nonunitary state, the product A(k)A(k) t = d(k) • d*(k) + *(d(k) x d*(k)) • a 
is not a multiple of unit matrix. Thus in a non-unitary state time reversal symmetry 
necessarily is broken spontaneously and spontaneous moment m(k) = id(k) x d*(k) 
appears at each point (k). In this case the macroscopic averaged moment < m(k) > 
integrated on the Fermi surface non-vanishes. The variable m(k) is related to the net 
spin average of state k by tr[A(k)"l"<3jA(k)], it is clear that the total spin average over 
the Fermi surface can be nonzero. As an application, the nonunitary bipolar state of 
/—wave pairing symmetry has been considered for the B— phase of superconductivity in 
the UPts compound which has been created in the low temperatures T and low values of 
H ^51 EI]- In this chapter, we want to investigate the weak link between two nonunitary 
superconducting bulks. This type of weak link structure can be used to demonstrate 
and to test the pairing symmetry in the superconducting phase (Mj. Consequently, 
we generalize the formalism of paper [E] for the case of weak link between nonunitary 
triplet superconducting bulks. In the paper [E|, the Josephson effect in the point contact 
between unitary /—wave triplet superconductors has been studied. Also, the effect of 
misorientation on the charge transport has been investigated and a spontaneous current 
tangential to the interface between the /—wave superconductors, has been observed. 
In this chapter, the ballistic Josephson weak-link via an interface between two bulks 
of nonunitary bipolar /—wave superconductivity with different orientations of the 
crystallographic axes is investigated. It is shown that the current-phase diagrams are 
totally different from the current-phase diagrams of the junction between conventional 
(s-wave) superconductors [30], high T c (d-wave) superconductors 48.] and unitary triplet 
( axial and planar) /—wave superconductors [E]. This different characters can be 
used to distinct between nonunitary bipolar /—wave superconductivity from the other 
types of superconductivity, roughly speaking. In this weak-link structure between the 
nonunitary /—wave superconductors, the spontaneous current parallel to the interface 
as a fingerprint of unconventional superconductivity and spontaneous time reversal 
symmetry breaking, has been observed. The effect of misorientation on the spontaneous 
and Josephson currents is investigated. It is possible to find the value of the phase 
difference in which the Josephson current is zero but the spontaneous current tangential 
to the interface, which is produced by the interface, is present. In some configurations 
and at the zero phase difference, the Josephson current is not zero generally and it has 
a finite value. This finite value corresponds to a spontaneous phase difference which is 
related to the misorientation between the gap vectors. 

The arrangement of the rest of this chapter is as follows. In Sec. (l6.2|) we describe our 
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configuration, which has been investigated. For a non-self-consistent model of the order 
parameter, the quasiclassiacl Eilenberger equations [3T] are solved and suitable Green 
functions have been obtained analytically. In Sec. 1)6.3)1 the obtained formulas for the 
Green functions have been used for calculation the current densities at the interface and 
an analysis of numerical results will be done. 



6.2. Quasiclassical equations 

We consider a model of a flat interface y = between two misoriented nonunitary 
/—wave superconducting half-spaces (FigH2J) as a ballistic Josephson junction. In 
order to calculate the current, we use quaisclassical Eilenberger equations [3T] for Green 
functions (J2J) The Matsubara propagator g here, like the unitary triplet superconductor 
case, can be written in the form: 



gx + gi • a g 2 • aia 2 
«0"2g3 • o ff 4 - o- 2 g 4 ■ oa 2 



(107) 



where, the matrix structure of the off-diagonal self energy A in the Nambu space is 

A=( ° . d '*M. (108) 
I ia 2 d* -a ) K ' 

In this chapter, the nonunitary states, for which dx d* ^ 0, is investigated. 
Fundamentally, the gap vector (order parameter) d has to be determined numerically 
from the self-consistency equation [S3], while in this chapter, we use a non-self-consistent 
model for the gap vector which is much more suitable for the analytical calculations. 
Solutions to Eq. (jl} must satisfy the conditions for Green functions and gap vector d 
in the bulks of the superconductors far from the interface as follow: 



where, 



and, 



1 




e(l-A a -a) 


id n - d n X A n 




( io-2 


[id* n + d* n x A n ] ■ o 


-£<7 2 (1 + A n - 






id n x d* n 






d n • d* + + d n ■ d*) 2 


+ (d n X d*) 2 


2 


(e 2 + d n ■ d n ) 2 + (d n 


X d*) 2 ] 



(109) 



V £ 2 + d n • d* + v/(£ 2 + d n • d*) 2 + (d n X d*) 2 

d(±oo) = d 2il (T, v F ) exp ( t4 ) (112) 



where <fi is the external phase difference between the order parameters of the bulks and 
n = 1, 2 label the left and right half spaces respectively. It is clear that, poles of Green 
function in the energy space, are in 

ih,! = (113) 



Charge and Spin Transport in Superconducting Weak Links 



60 




Figure 32. Component of current normal to the interface (Josephson current) versus 
the phase difference <fi for the junction between nonunitary bipolar /—wave bulks , 
T/T c = 0.15, geometry (i) and different misorientations. Currents are given in units 
of jo = fe7V(0)« F A (0). 



consequently, 

((zEf + d 2)1 ■ d^) 2 + (d 2)1 x d* 2il ) 2 = (114) 

and 

E = ±y / d 2il -d^ lT zd 2>1 xd^ 1 (115) 

in which E is the energy place of poles. Eq. have to be supplemented by the 
continuity conditions at the interface between superconductors. For all quasiparticle 
trajectories, the Green functions satisfy the boundary conditions both in the right and 
left bulks as well as at the interface. 

The system of equations (0) and the self-consistency equation |HS] can be solved only 
numerically. For unconventional superconductors such solution requires the information 
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of the interaction between the electrons in the Cooper pairs and nature of unconventional 
superconductivity in novel compounds which, in most cases are unknown. Also, It has 
been shown that the absolute value of a self-consistent order parameter is suppressed 
near the interface and at the distances of the order of the coherence length, while its 
dependence on the direction in the momentum space almost remains unaltered |41j . 
This suppression of the order parameter which keeps the current-phase dependence 
unchanged but changes the amplitude value of the current, doesn't influence the 
Josephson effect drastically. For example, it has been verified in paper [48] for the 
junction between unconventional cf-wave, in paper jH] for the case of unitary "/- 
wave" superconductors and in paper [Slj for pinholes in 3 He that, there is a good 
qualitative agreement between self-consistent and non-self-consistent results. Also, it 
has been observed that results of the non-self-consistent model in are similar to the 
experiment [13] • Consequently, despite the fact that this estimation cannot be applied 
directly for a quantitative analyze of the real experiment, only a qualitative comparison 
of calculated and experimental current-phase relations is possible. In our calculations, 
a simple model of the constant order parameter up to the interface is considered and 
the pair breaking and the scattering on the interface are ignored. We believe that under 
these strong assumptions our results describe the real situation qualitatively. In the 
framework of such model, the analytical expressions for the current can be obtained for 
a certain form of the order parameter. 

6.3. Analytical and numerical results 

The solution of Eq. allows us to calculate the current densities. The expression for 
current is: 



where, (...) stands for averaging over the directions of an electron momentum on the 
Fermi surface v> and N (0) is the electron density of states at the Fermi level of energy. 
We assume that the order parameter is constant in space and in each half-space it equals 
to its value (|112j) far form the interface in the left or right bulks. For such a model, 
the current-phase dependence of a Josephson junction can be calculated analytically. 
It enables us to analyze the main features of current-phase dependence for a certain 
model of the order parameter of nonunitary /—wave superconductivity (bipolar). The 
Eilenberger equations (pQ) for Green functions g, which are supplemented by the condition 
of continuity of solutions across the interface, y = 0, and the boundary conditions at 
the bulks, are solved for a non-self-consistent model of the order parameter analytically. 
In a ballistic case the system of equations for functions gi and gj can be decomposed 
on independent blocks of equations. The set of equations which enables us to find the 
Green function g\ is: 




(116) 



in 



v F k\7g! = i (d • g 3 - d* • g 2 ) ; 
v F kVg- = -2 (d x g 3 + d*xg 2 ) ; 



(117) 
(118) 
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Figure 33. Component of current normal to the interface (Josephson current) versus 
the phase difference <fi for the junction between nonunitary bipolar /—wave bulks , 
T/T c = 0.15, geometry (ii) and different misorientations. 



u F kVg 2 = -2e m g 2 + 2igid + d x g_; (119) 
v F kVg 3 = 2e m g 3 - 2i 9l d* + d* x g_ ; (120) 

where g_ = gi — g 4 . The Eqs. (jll7|) - (|119|) can be solved by integrating over ballistic 
trajectories of electrons in the right and left half-spaces. The general solution satisfying 
the boundary conditions ()109|) at infinity is 

Si n) = ^ + ane- 2sfW ; (121) 
g W = -2^A n + C n e- 2s ^; (122) 

g2 " 2sr ] n n -2e m 6 ] (123) 
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(n) 



+ d* 

1 n 



x A 



n 



+ 



2ia n d* n - d; x 
2sr]Vt n + 2s. 



x C 



2 S^&yz ~t 



(124) 



rn 



where t is time of flight along the trajectory, sgn (t) = sgn (y) = s and r] = sgn (v y ) . 
By matching the solutions (jl21H124j) at the interface (y = 0, t = 0), we find constants a n 
and C n . Indices n = 1,2 label the left and right half-spaces respectively The function 
9i (0) = g^p (— 0) = g^ (+0) , which is a diagonal term of Green matrix and determines 
the current density at the interface, y = 0, is as follows: 



in the momentum space. The crystallographic c-axis in the left half-space is selected 
parallel to the partition between the superconductors (along z-axis in FiglT^j). To 
illustrate the results obtained by computing the formula (]125j) . we plot the current- 
phase diagrams for two different geometries. These geometries are corresponding to the 
different orientations of the crystals in the right and left sides of the interface (Fig^J): 

(i) The basal afe-plane in the right side has been rotated around the c-axis by a; Ci||c 2 . 

(ii) The c-axis in the right side has been rotated around the 6-axis by a; bi||b 2 . 
Further calculations require a certain model of the gap vector (order parameter) d. 

In this chapter, the nonunitary /—wave gap vector in the B— phase (low temperature 
T and low field H) of superconductivity in UPt^ compound has been considered. This 
nonunitary bipolar state which explains the weak spin-orbit coupling in UPt^ is [15] : 



left side of FigH21 The function A =A (T) describes the dependence of the gap 
vector on the temperature T (our numerical calculations are done at the low value of 
temperature T/T c = 0.1). Using this model of the order parameter (|126p and solution to 
the Eilenberger equations ()125]) . we have calculated the current density at the interface 
numerically. These numerical results are listed below [57] : 

1) The nonunitary property of Green's matrix diagonal term consists of two part. The 
explicit part which is in the B mathematical expression in Eq. ([125]) and the implicit part 
in the fii^ and d 12 terms. These f2i i2 and d 12 terms are different from their unitary 
counterparts. In the mathematical expression for f2i j2 the nonunitary mathematical 
terms A 12 are presented. The explicit part will be present only in the presence of 
misorientation between gap vectors, 

B = idi x d 2 • (Ax+A 2 )(?7fi 2 — e)(r]Q 1 +s), but the implicit part will be present always. 
So, in the absence of misorientation (di||d 2 ), although the implicit part of nonunitary 
exists but the explicit part is absent. This means that, in the absence of misorientation 
current-phase diagrams for planar unitary and nonunitary bipolar systems are the same 
but the maximum values is different slightly. 

2) A component of current parallel to the interface j z for geometry (i) is zero as same as 




The coordinate axes x 
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0.25 0.5 0.75 1 

Figure 34. The x— component of charge current tangential to the interface 
versus the phase difference </> for the junction between nonunitary bipolar /—wave 
superconducting bulks, T/T c = 0.15, geometry (i) and the different misorientations. 

the unitary case jE] while the other parallel component j x has a finite value (see FiglMj). 
This later case is a difference between unitary and nonunitary cases. Because in the 
junction between unitary /—wave superconducting bulks all parallel components of the 
current (j x and j z ) for geometry (i) are absent 0. 

3) In FigsE2ESl the Josephson current j y is plotted for certain nonunitary model of 
/—wave and different geometries. Figs l32l33l are plotted for the geometry (i) and 
geometry (ii) respectively. They are completely unusual and totally different from their 
unitary counterparts which have been obtained in j^]. 

4) In Fig|H21for geometry (i), it is observed that by increasing the misorientation, some 
small oscillations, as the result of non-unitary property of the order parameter, appear 
in the current-phase diagrams. Also, the Josphson current at the zero external phase 
difference <fi = is not zero but it has a finite value. The Josephson current will be zero 
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Figure 35. Current tangential to the interface versus the phase difference <fi for the 
junction between nonunitary bipolar /—wave superconducting bulks, T/T c = 0.15, 
geometry (ii) and the different misoricntations (^component) . 



at the some finite values of the phase difference. 

5) In FigEmfor geometry (ii), it is observed that by increasing the misorientation the 
new zeros in current-phase diagrams appear and the maximum value of current will 
be changed non-monotonically. In spite of the FidH21 for geometry (i), the Josephson 
current at the phase differences = 0, = tt and <fi = 2ir is zero exactly. 

6) The current-phase diagram for geometry (i) and x— component fFiglMj) is totally 
unusual. By increasing the misorientation, the maximum value of current increases. 
The components of current parallel to the interface for geometry (ii) are plotted in 
Fig Ellwand FielSBI All of terms at the phase differences <fi = 0, <fi = tt and <fi = 2tt are 
zero. The maximum value of the current-phase diagrams is not a monotonic function of 
the misorientation. 
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Figure 36. Current tangential to the interface versus the phase difference <fi for the 
junction between nonunitary bipolar /—wave superconducting bulks, T/T c = 0.15, 
geometry (ii) and the different misorientations (z— component). 



7. Conclusions 



At first, in the Chapter 3 of this thesis, we have studied the stationary Josephson effect 
in the ballistic point contact with transport current on the banks in the model S — c — S 
taking into account the reflection of electrons from the contact. The contact is subject 
to the phase difference <fi and the transport current tangential to the boundary of the 
contact. As it was shown in in the ideal transparent point-contact at <ft = it and 
near the orifice the tangential current flows in the opposite direction to the transport 
current, and there are two anti-symmetric vortex-like structures. It is observed that, 
by decreasing the transparency D c < D < 1 the vortex-like current is destroyed. The 
critical value of D = D C (T) depends on the temperature T and D c (T — > 0) — > 1, 
D c (T — > T c ) — > |, so that we can never find a vortex for transparency values lower than 
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| . This anomalous temperature behavior of the vortices is the result of non-monotonic 
dependence of the interference current on the temperature. 

Then, we have studied the spin current in the ballistic Josephson junction in the model of 
an ideal transparent interface between two misorientated /-wave superconductors which 
are subject to a phase difference 0, in Chapter (@J). Our analysis has shown that the 
misorientation and different models of triplet gap vectors influence the spin current [83J . 
This has been shown for the charge current in the paper Misorientation between 
the left and right gap vectors changes strongly the critical values of both of the spin 
current and charge current. It has been obtained that the spin current is the result of the 
misorientation between the gap vectors. Furthermore, as an interesting and new result, 
it is observed that the different models of the gap vectors and geometries can be applied 
for the polarization of the spin transport. Another result of these calculations is the state 
in which the currents select one of the two possible directions (perpendicular and parallel 
to the interface) to flow. This property can be used as a switch to control the direction 
of the charge and spin current. Finally, it is observed that at some certain values of the 
phase difference 0, the charge-current vanishes while the spin-current flows, although 
the carriers of both spin and charge are the same (electrons). The spatial variation 
of the phase of the order parameter plays a role as the origin of the charge current 
and, similarly, due to a spatial difference of the gap vectors in two half-spaces causes 
spin currents. This is because there is a position-dependent phase difference between 
"spin up" and "spin down" Cooper pairs and, although the total charge current vanishes, 
there can be a net transfer of the spin. Therefore, in our system, there is a discontinuous 
jump between the gap vectors and, consequently the spin currents should generally be 
present. For instance, if spin-up states and spin-down states have a velocity in the 
opposite directions, the charge currents cancel each other whereas the spin current is 
being transported. Mathematically speaking, j c har g e = j| + j|Js P in = jt — jj. , so it is 
possible to find the state in which one of these current terms is zero and the other term 
has a finite value jHBl- m addition, the spin imbalance which is the result of the different 
density of states for "spin-up" and "spin-down" can be other reason of spin current jHI|. 
In conclusion, the spin current in the absence of the charge current can be observed. 
Also, the weak link between two misorientated PrOs^Sbu superconducting bulks has 
been investigated in Chapter (JHJ) of this thesis. We have considered the transportation of 
the spin and charge in the ballistic Josephson junction between two PrOs^Sb^ crystals 
with u (p+h)— wave" pairing symmetry which are subject to a phase difference 0. In this 
case as same as the case of /—wave junctions in Chapter the different misorientations 
and different models of the gap vectors influence the spin and charge currents and it 
is shown that the misorientation of the superconductors leads to a spontaneous phase 
difference that corresponds to the zero Josephson current. This phase difference depends 
on the misorientation angle. A spontaneous charge current tangential to the interface 
which is not equal to zero in the absence of the Josephson current is observed in this 
junction. It has been obtained that the spin current is the result of the misorientation 
between the gap vectors which is the characteristic of unitary triplet superconductors. 
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Again in this case as same as the case of /—wave superconductors, it is observed that 
the certain model of the gap vectors and geometries can be applied to polarize of the 
spin transport. Finally, it is observed that at some certain values of the phase difference 
0, the charge-current vanishes while the spin-current flows. The reason for this case 
is the above-mentioned discussion for /—wave. The spin and charge currents can be 
used to recognize the A— phase and B— phase of "(p + h)— wave" and the pure p— wave 
pairing symmetry and can be used to determine the pairing symmetry. Particularly, this 
proposed experiment can be used to demonstrate the h)— wave" pairing symmetry 
for which many doubts and challenges exist. 

Finally, in Chapter (jZ|), we have theoretically studied the charge currents in the 
ballistic Josephson junction in the model of an ideal transparent interface between two 
misoriented U Pt$ crystals with nonunitary bipolar /—wave superconducting bulks which 
are subject to a phase difference <fi. Our analysis has shown that misorientation between 
the gap vectors create a current parallel to the interface and the different misorientations 
between gap vectors influence the spontaneous parallel and normal Josephson currents. 
These have been shown for the currents in the point contact between two bulks of 
unitary axial and planar /-wave superconductor in jH] separately. Also, It is shown that 
the misorientation of the superconductors leads to a spontaneous phase difference that 
corresponds to the zero Josephson current and to the minimum of the weak link energy 
in the presence of the finite spontaneous current. This phase difference depends on the 
misorientation angle. Again in this case, the tangential spontaneous charge current is not 
equal to zero in the absence of the Josephson current. The difference between junction 
behavior of unitary planar and nonunitary bipolar superconductivity can be used to 
distinct between them. This experiment can be used to test the pairing symmetry and 
recognize the different phases of UPt^ which has two unitary and a nonunitary triplet 
pairing symmetry. 
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